NEW  MEXICO  STATE  UNIV  LAS  CRUCES  PHYSICAL  SCIENCE  LAP  F/G  4/1 

D-REGION  ELECTRON  DENSITY  MEASUREMENTS  OBTAINED  BY  PARTIAL  REFL— ETC(U) 
NOV  60  D  MOTT#  W  GAMMILL#  R  VALDEZ  DAAD07-79-C-0008 

ERADCOM/ASL-CR-0O-OOO8-4 


.1 


w 


DOG  FILE  COP^  ADA095566 


D-REGION  ELECTRON  DENSITY 
MEASUREMENTS  OBTAINED  BY 
PARTIAL  REFLECTION  TECHNIQUES 

NOVEMBER  1980 


Prepared  by 

l>.  MOTT 
W.  GAMMILL 
R.  VALDEZ 

Physical  Science  Laboratory 
New  Mexico  State  University 
Las  Cruces,  New  Mexico  88003 


UNDER  CONTRACT:  DAAD07-79-C0008 

Contract  Monitor:  Robert  Olsen 


Approved  for  public  release;  distribution  unlimited 


US  Army  Electronics  Research  and  Development  Command 

ATMOSPHERIC  SCIENCES  LABORATORY 

White  Sands  Missile  Range,  NM  88002 


81  2  26  047 


NOTICES 


Di scl aimers 


The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig¬ 
nated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  be  construed  as  official  Government 
indorsement  or  approval  of  commercial  products  or  services 
referenced  herein. 


Disposition 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not 
return  it  to  the  originator. 


S  ECU  RITYCLASSI  FI  CATION  OF  THIS  PAGE  (Whan  Data  Entered) 


C/9)  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

ASL-fCMf  0008-4  ]  ^ f09S~  Si 

3.  RECIPIENT’S  CATALOG  NUMBER 

(9); 

4.  TITLE  Cmd  SubllU*) _ — 

J1-.BEGI0NJLECTR0NJJENSITY  MEASUREMENTS  *- 
OBTAINED  BY  PARTIAL  REFLECTION  TECHNIQUES#  l 

5.  TYpI^O’THEPORT  ft  PEWODTCOVERED 

Mtif Final  Re^erJ^  j 

_  *  \ 

|  -  ’•  ]  _ 

-r-nrrymnfm)  j  ,  ,  ( 

D.  Mott/ W.y6ammil  1 /R./Valdez  j 

ft.  CONTRACT  OR  GRANT  NUMBERfmJ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Physical  Sciences  Laboratory 

New  Mexico  State  University 

Las  Cruces,  New  Mexico  88003 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA^HMORK  UNIT  NUMBERS 

—  7 

DA  Ta  sic/  DAAD07-79-9p008 / ^ 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  - V 

US  Army  Electronics  Research  (//) 

12-  REEOJtT-D  AT  E 

Novamkmr  W80/ 

and  Development  Command 

Adel phi,  MD  10783 

T«.  MONITORING  AGENCY  NAME  ft  AOORESSff/  different  from  Controlling  Ottlca) 

Atmospheric  Sciences  Laboratory 

White  Sands  Missile  Range,  IN  88002 

is.  NUMBER  OF  PAGES 

125 

IS.  SECURITY  CLASS,  (ol  thtm  nport) 

UNCLASSIFIED 

15*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (ol  (A It  Roport) 

Approved  for  public  release;  distribution  unlimited 

t  •'  j  ^  /  ✓ 

.  z 

*•-  SUPPLEMENTARY  NOTES 

Robert  Olsen,  Contract  Monitor 

It*  KEY  WOROS  (Continue  on  taeara •  aid*  it  nacMMry  and  identify  by  block  number) 

D- region  HE  propagation 

Electron  densities  Refractive  Index 

Partial  reflection  sounder  Ionosphere 


|Mt  ABSTRACT  CGhBm  am  tavarwa  alM  «#  name  meaty  mtd  tdentify  by  block  numbar) 

^D-reglon  electron  densities  are  derived  utilizing  a  ground  based  partial 
reflection  sounder.  The  sounder  frequency  Is  approximately  2  MHz,  with 
right-hand  and  left-hand  polarized  pulses  transmitted  vertically.  This  report 
describes  the  sounder  system  hardware  and  the  necessary  software  to  derive 
electron  densities  from  the  returned  vertical  pulses.  A  complete  set  of 


DO  MTn  un  EDITION  OP  t  NOV  SS  IS  OMOLETC 


security  classification  of  this  pace  r« 


A' 

i  Data  Bn  tor  ad) 


20.  ABSTRACT  (cont) 


schematics  of  the  electronics,  antenna  array  and  operational  are  Included.  A 
total  of  604  electron  density  profiles  have  been  reduced  with  the  sounder 
located  at  White  Sands  Missile  Range  (WSMR)  and  at  Red  Lake,  Canada,  as  part 
of  the  1979  solar  eclipse  field  program.  The  WSMR  profiles  have  been 
summarized  Into  monthly  profiles.  The  electron  densities  derived  from 


utilizing  the  partial  reflection  technique  appear  to  be  within  the  expected 
ranges  for  D-reglon  profiles. 


ttCUftlTY  CLASSIFICATION  OF  THIS  PAOCOWimfl  Dmtm  Knfrmd) 


CONTENTS 


SECTION 

1.  Theory  of  the  partial  reflection  experiment .  7 

2.  Experimental  procedure  .  12 

3.  System  description  .  14 

4.  Antenna  installation  .  25 

5.  System  operation  .  29 

6.  Data  processing . 35 

7.  An  alternate  procedure  for  data  processing . 49 

8.  Results . 53 

REFERENCES . 76 

APPENDIX 

I.  Circuit  diagrams  77 

II.  Data  averaging  computer  program  r  107 

III.  Solar  Zenith  angle  115 

IV.  Receiver  calibration  120 

V.  Polynomial  approximations  for  auc*  C5/2(x)  124 


3 


LIST  OF  FIGURES 


3.1  Partial  reflection  sounder,  block  diagram  .  15 

3.2A  Exterior  photograph  of  van .  16 

3.2B  Interior  photographs  of  van .  17 

3.3  Phases  for  transmitting  with  circular  polarization .  19 

3.4  Phases  for  receiving  with  circular  polarization  .  21 

4.1  Dipole  antenna  element .  25 

4.2  Antenna  array,  plan  view .  26 

4.3  Antenna  matching  calibration  setup .  28 

5.1  Receiver  calibration  setup .  33 

*6.1  Accumulate  a  logic  for  data  averaging .  36 

6.2  Output  of  averaging  program  .  39 

6.3  Procedures  used  in  standard  data  processing .  44 

6.4  Example  result,  standard  data  processing .  48 

7.1  Example  result,  non-standard  data  processing .  52 

8.1  Key  to  electron  density  profiles .  54 

8.2  WSMR  data  averages,  September  1977 .  65 

8.3  WSMR  data  averages,  October  1977 .  66 

8.4  WSMR  data  averages,  November  1977  .  67 

8.5  WSMR  data  averages,  December  1977  .  68 

8.6  WSMR  data  averages,  January  1978 .  69 

8.7  WSMR  data  averages,  February  1978  .  70 

8.8  WSMR  data  averages,  March  1978 .  71 

8.9  WSMR  data  averages,  April  >978 .  72 

8.10  WSMR  data  averages,  August  1978  .  73 

8.11  WSMR  data  averages,  July  1978  .  74 

8.12  Electron  density  vs.  time;  solar  eclipse  of  February  1979  ...  75 

AIII-1  Solar  zenith  angle  triangle  .  115 

AIII-2  Solar  zenith  angle  computer  program  .  116 


4 


LIST  OF  TABLES 


5.1  Computer  bootstrap  procedure . 31 

6.1  Amplitude  calculations . 42 

6.2  Tabulated  functions  R(h)  and  G(h) . 45 

6.3  Specifications  for  Table  6.2 . 46 

7.1  Tabulated  functions  E(h)  and  u(h) . 51 

8.1  Run  index  for  WSMR,  1975 .  56 

8.2  Run  index  for  WSMR*,  1977-1979  .  57 

8.3  Run  index  for  Balmertwon,  Canada,  1979 .  63 

AIII-1  Base  values  for  zenith  angle  calculations  .  119 

AIV-1  Typical  receiver  calibration  data  .  121 

AIV-2  Tabulated  receiver  calibration . 123 


5 


SECTION  1 


Theory  of  the  Partial  Reflection  Experiment 

This  experiment  is  designed  to  measure  the  density  of  free  electrons  in  the 
ionosphere  as  a  function  of  altitude,  in  the  D  region  below  100  kilometers.* 

A  radar  operating  at  a  frequency  of  several  megahertz  transmits  short  pulses 
of  radiation  vertically,  and  echoes  are  partially  reflected  (backscattered) 
from  the  free  electrons  in  the  lower  ionosphere.  Echo  strength  vs.  time  (with 
time  a  measure  of  altitude)  is  digitized  and  recorded  on  magnetic  tape. 

Circular  polarization  of  the  electromagnetic  radiation,  both  right-hand  and 
left-hand  polarization,  is  utilized.  Several  pulses  are  transmitted  and 
received  in  right-hand  polarization,  and  the  next  several  pulses  are  trans¬ 
mitted  and  received  in  left-hand  polarization. 

Because  of  the  earth’s  magnetic  field,  the  index  of  refraction  of  the  ionosphere 
is  different  for  the  two  polarization  modes.  This  difference  leads  to  different 
coefficients  of  reflection  and  absorption  of  the  radio  waves.  The  mode  for 
which  absorption  and  reflection  is  less  is  called  the  ordinary  mode,  while  the 
other  is  called  the  extraordinary  mode.  The  differences  in  the  absorption  and 
reflection  provide  the  key  to  measuring  electron  density. 

It  is  assumed  that  the  amplitude  Am  (m  is  an  index  to  identify  polarization 
mode)  of  an  echo  partially  reflected  from  the  ionosphere,  from  altitude  h 
above  the  station,  is  related  to  an  ionospheric  reflection  coefficient  Rm  and 
absorption  coefficient  K  as  follows: 


( k)  oc  Rfm(h)  e 


As  shown  below,  Km  at  any  altitude  is  proportional  to  the  density  N  of  free 

electrons  at  that  altitude;  i.e.,  K  =  F  N,  where  F  is  the  proportionality 

mm  m 

factor . 


*It  is  also  known  as  the  differential  absorption  experiment.  Reference  (1) 
contains  an  excellent  review  of  the  theory. 
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When  circularly  polarized  waves  of  both  the  ordinary  (o)  mode  and  the  extraor¬ 
dinary  (x)  mode  are  transmitted  at  equal  power,  then 


t. 


This  equation  is  solved  for  N: 


N00 


i 

^[Fx-Folv, 


Ro 

Ak/ a0 


(1.1) 


Equation  (1.1)  may  serve  as  the  basis  for  reducing  data  from  the  partial  re¬ 
flection  experiment.  In  this  case,  the  function  A^/Ao  vs.  altitude  is  the 
input  data,  while  functions  R^/R^  and  [F^  -  Fq]  are  precomputed  from  theo¬ 
retical  relationships.  The  reader  should  note  that,  in  Equation  (1.1),  the 
mathematical  operations  performed  on  the  data  are  such  that,  under  some  condi¬ 
tions,  a  small  fractional  uncertainty  in  or  Aq  may  result  in  a  large  frac¬ 
tional  uncertainty  in  N. 

According  to  Sen  and  Wyller  (2),  and  in  the  longitudinal  approximation  (direction 
of  propagation  parallel  to  the  earth’s  magnet  field  lines),  the  square  of  the 
ionosphere's  index  of  refraction,  as  a  function  of  altitude,  is  this: 

cl(  yJ  +iic^*)]  (1,2) 


oo  V 


’+'  corresponding  to  ordinary  circular  polarization  of  the  electromagnetic 
radiation,  and  to  the  extraordinary  polarization.  Equation  (1.2)  contains 
these  terms: 
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Me 

N(K) 
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plasma  frequency,  radians/sec 
density  of  free  electrons  at  altitude  h 
electron  charge 
electron  mass 


£0 

60 

=  eS 
"  /m. 

B 


permittivity  of  free  space 

signal  frequency,  radians/sec 

gyromagnetic  frequency,  radians/sec  (assumed 
constant  over  altitudes  of  interest) 

earth* s  magnetic  field  intensity 


iXh) 


CfM  a 


f  -t 


V  e 
/x*+  ex 


-4  £ 


rate,  in  collisions  per  second,  at  which  free 
electrons  at  altitude  h,  having  the  most  probable 
energy,  collide  with  other  particles 

See  Appendix  V  for  numerical  approximations  of 
this  integral. 


It  is  well  known  that  the  reflection  coefficient  for  the  ionosphere  is 
proportional  to  the  fractional  variation  in  n,  i.e.,  R  oc  b/n/n  .  There  is 
general  agreement  among  investigators  that  variations  in  n  are  proportional 
to  variations  in  electron  density,  and  that  n  is  relatively  independent  of 
variations  in  other  parameters  of  the  ionosphere.  Accordingly,  S/>tacSN  is 
assumed.  From  Equation  (1.2)  it  then  follows  that 


R  CC - - — 

'm  /n*  \ 


In  the  partial  reflection  experiment,  only  the  absolute  value  of  the  reflection 
coefficient  has  physical  significance.  And  in  the  D  region  of  the  ionosphere, 
n+  and  n  both  have  approximately  unit  value.  Thus,  the  function  R^/R^  in 
Equation  (1.1)  may  be  evaluated  from  this  expression: 


ROO  s 


L y.  f-)]  +  tz 


(1.3) 
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The  ionosphere's  index  of  refraction  has  a  real  part  {jm  )  and  an  imaginary 
part  (?£)>  i.e.,  .  The  imaginary  part  is  proportional  to  the 

absorption  coefficient:  K  =  X  .  Now  /n2  =  - i  7.^ X  ,  and  so  by 

comparing  this  expression  with  Equation  (1*2)  we  see  that 

K~=  If  =  M  , 

where  JJLm*  f  has  been  assumed.  Thus,  the  function  2(Fx-Fq)  in  Equation  (1.1) 
may  be  evaluated  from  this  equation: 

GOn)=  2(^-0  =  T  C|^t+)j  (j  4) 


Flood  (3)  has  developed  an  improved  approximation  for  index  of  refraction  that 
is  comparable  to  the  Sen  Wyller  longitudinal  approximation  of  Equation  (1.2), 
but  that  is  more  accurate  when  the  propagation  is  not  strictly  longitudinal. 
His  expression  corresponding  to  Equation  (1.2)  is  this: 


/n 


z 

± 


(*>H 

CO  V 


A  L  ^±Ci(^±)  +  i  x  C|.(^±)  ] 


+  B[«*C*Ofc>  +  i  4  <*(•**)  ] 


( !■  2  A) 


whare 


u> 

it=  — 

A  =  cos*  4*  ~  ^  s,r,a^ 
B  =  -  -jf  Sin1^ 

i>  =  ~k 


10 


with  <f>  the  angle  between  the  magnetic  field  lines  and  the  direction  of  propa¬ 
gation  (which  is  vertical  in  our  case).  All  other  terms  here  are  the  same  as 
in  Equation  (1.2).  Note  that  Equation  (1.2A)  reduces  to  Equation  (1.2)  in  the 
special  case  (f  —  0. 


One  progresses  from  Equation  (1.2A)  to  the  corresponding  new  functions  for 
R(h)  and  G(h)  using  the  same  logic  as  that  leading  from  Equation  (1.2)  to 
Equations  (1.3)  and  (1.4).  Results  are  as  follows: 


R(W)  = 


]  +  (a)  1 

1 I A  %  ci ( B f  ci  ^-)+D9ci  (*)] 

)]2  J 

x 


icy 


(1.3A) 

(1.4A) 


In  summary,  the  reduction  of  partial-reflection  data  utilizes  Equation  (1.1)  to 
convert  the  ratios  (vs.  altitude)  of  the  measured  amplitudes  of  extraordinary 
and  ordinary  echoes,  A^/A^,  into  an  electron  density  vs.  altitude  profile. 
Equation  (1.1)  requires,  in  addition  to  the  amplitude  ratios,  precomputed 
functions  R(h)  and  G(h)  from  Equations  (1.3A)  and  (1.4A).  These  functions  are 
evaluated  using  either  an  experimental  or  theoretical  profile  of  collision 
frequency  v(h)  ,  and  an  experimental  or  theoretical  value  of  the  gyromagnetic 
frequency  oUH  at  D-region  altitude.  Using  the  R(h)  and  G(h)  notation  from 
Equations  (1.3)  and  (1.4),  Equation  (1.1)  may  be  written  in  this  form: 


d  h 


(1.5) 


Tabulations  of  R(h)  and  G(h)  used  in  this  study  are  given  in  Table  6.2. 
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SECTION  2 


Experimental  Procedure 

In  practice,  the  amplitudes  and  Aq  (vs.  altitude)  of  the  extraordinary  and 
ordinary  echoes  that  are  used  in  Equation  (1.5)  are  average  values.  Echo 
amplitudes  are  subject  to  rapid,  random  fluctuation  (4),  and  if  one  were  to 
use  in  Equation  (1.5)  echo  amplitudes  from  single  measurements,  the  result 
would  not  be  statistically  significant.  The  time  of  a  typical  "run"  of  the 
experiment  is  10  minutes,  and  the  averages  extend  over  such  a  time  interval. 
Electron  density  in  the  lower  ionosphere  is  generally  a  function  of  solar 
zenith  angle,  and  the  change  in  this  angle  over  a  10-minute  period  is  negli¬ 
gible  . 

The  partial-reflection  transmitter  operates  at  a  pulse  repetition  frequency  of 
about  17  pulses  per  second,  and  so  a  great  deal  of  data  is  accumulated  in  the 
course  of  a  run.  All  of  this  data  is  recorded  on  magnetic  tape  for  subsequent 
analysis . 

Eight  consecutive  pulses  from  the  transmitter  form  the  basi<~  pattern  of  system 
operation,  which  is  under  computer  control.  During  the  first  4  pulses,  trans¬ 
mission  and  reception  is  done  in  the  ordinary  mode,  and  during  the  second  4 
pulses,  in  the  extraordinary  mode.  During  each  group  of  4  pulses,  receiver 
attenuation  is  increased  in  steps,  with  no  attenuation  for  the  first  pulse, 
one  unit  of  attenuation  for  the  second  pulse,  two  units  for  the  third,  and 
three  units  for  the  fourth.  The  operator  selects  the  unit  of  attenuation 
here,  which  may  be  either  6,  9,  or  12  dB. 

Four  levels  of  receiver  attenuation  are  used  to  insure  that  some  echoes  from 
the  altitude  region  of  interest  are  in  the  range  of  receiver  sensitivity; 
above  the  noise  level  but  below  the  saturation  level.  In  addition  to  the  4 
stepped  levels  of  attenuation,  the  operator  can  add  fixed  attenuation  at  the 
receiver  input,  and  this  increases  the  attenuation  at  all  levels  by  the  same 
amount.  As  an  example,  suppose  that  the  unit  of  attenuation  for  the  stepping 
is  selected  to  be  6  dB.  Then  the  four  levels  of  stepped  attenuation  are  0, 
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6  dB,  12  dB,  and  18  dB.  But  if  the  operator  selects  a  fixed  value  of  attenua¬ 
tion  of  10  dB,  then  the  net  attenuation  in  the  4  steps  is  10  dB,  16  dB,  22  dB, 

and  28  dB. 

The  echo  from  each  pulse  is  "stretched"  because  of  the  partial  reflections 
that  occur  from  all  altitudes  in  the  lower  ionosphere.  That  is,  the  extent  in 

time  of  the  echo  is  much  greater  than  that  of  the  transmitted  pulse,  and,  in 

fact,  the  shape  of  the  echo  pulse  is  the  information  of  interest.  Each  echo 
is  digitized  for  recording  on  magnetic  tape  so  that  the  amplitude  functions 
A^Ch)  and  AQ(h)  may  be  determined.  The  current  digitizer  is  a  6-bit  (only) 
device,  so  each  value  is  an  integer  between  0  and  63. 

The  transmitted  pulse  has  a  time  duration  of  20  x  10  6  sec,  making  its  spatial 
extent  6  kilometers  (pulse  time  multiplied  by  the  speed  of  light).  It  follows 
that  the  echoes  being  received  at  one  instant  may  have  been  reflected  from 
anywhere  within  a  zone  3  kilometers  thick,  and  variations  in  electron  density 
with  altitude  that  occur  on  a  scale  smaller  than  this  will  not  be  detected. 
Accordingly,  each  echo  is  digitized  at  a  sample  rate  corresponding  to  a  sepa¬ 
ration  distance  Ah  of  2  km  between  samples,  and  30  samples  per  echo  are  recorded. 
Thus  echo  information  from  an  altitude  interval  of  60  km  is  recorded.  The 
base  of  this  altitude  interval  is  selected  by  the  operator  so  that  the  "leading 
edge”  of  the  echoes  is  always  sampled.  The  lower  altitude  is  set  typically  to 
between  50  and  60  kilometers.  If  50  km  is  chosen,  then  the  first  echo  sample 
corresponds  to  return  (if  any)  from  50  km  (with  altitude  uncertainty  determined 
by  pulse  width),  and  the  30th  and  last  sample  corresponds  to  return  from  110 
km.  Results  from  the  partial  reflection  experiment  are  generally  not  reliable 
above  about  90  km  due  to  the  fact  that  most  of  the  transmitted  energy  tends  to 
be  absorbed  or  reflected  at  lower  altitudes,  especially  in  the  extraordinary 
mode . 
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SECTION  3 


System  Description 

In  this  section,  the  components  of  the  partial  reflection  sounder  are  described. 
The  block  diagram  of  the  system  is  given  in  Figure  3.1.  Electronic  equipment 
is  mounted  in  a  mobile  van,  views  of  which  are  shown  in  Figure  3.2. 

3.1.  Power  Supply  and  Relay  Drivers 

This  unit  supplies  ±5  Vdc  and  ±15  Vdc  to  the  control  logic.  The  -5  V  level 
provides  a  negative-logic  capability  not  presently  used.  The  power  supply 
also  provides  -24  Vdc  to  the  relay  drivers  that  control  antenna  polarization. 

3.2.  Control  Logic 

The  control  logic  unit  is  the  source  of  all  commands  to  the  transmitter  and 
receiver . 

The  following  parameters  are  selected  by  the  operator  using  the  control  unit: 


Parameter 


Normal  Value 


pulse  repetition  frequency, 
17,12,  or  8  pps 


17  pps 


pulse  width 

sample  starting  height 

attenuation  stepping  unit, 
6,  9,  or  12  dB 


20  microseconds 
between  50  and  60  km 
6  dB 


Antenna  polarization 


0  and  X  mode 


Switches  on  the  control-unit  panel  are  used  to  enter  status  information  (date, 
time,  etc.)  that  is  recorded,  along  with  data,  on  the  magnetic  tape. 
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Figure  3.  !•  Partial  reflection  sounder  block  diagram 


3.3.  Pulse  Generator 


The  pulse  generator  is  used,  under  normal  conditions,  only  to  amplify  the  prf 
waveform  from  the  control  logic.  Normal  output  level  is  approximately  40 
volts  peak  to  peak,  which  is  sufficient  to  drive  the  class  C  exciter.  The 
generator  is  manufactured  by  E-H  Research  Laboratory,  model  131. 

3.4.  Transmitter 

The  transmitter  is  divided  into  three  parts;  exciter,  driver,  and  power  ampli¬ 
fier.  Each  part  operates  with  class  C  bias.  Four  switches  control  transmitter 
power,  as  follows: 

51  -  filament  voltage  to  exciter  and  driver; 

52  -  high  voltage  to  exciter; 

53  -  filament  voltage  to  power  amplifier  (P.A.); 

54  -  high  voltage  to  driver  and  P.A. 

Normal  operating  plate  voltage  on  the  P.A.  is  12.5  kV.  A  variac  is  used  to 
set  this  voltage  level.  With  12.5  kV  on  the  plates,  the  P.A.  screens  receive 
2.7  kV. 

There  is  a  90-sec  electronically-controlled  delay  after  P.A.  filaments  are 
turned  on,  before  the  P.A.  high  voltage  can  be  turned  on.  It  is  advisable  to 
let  the  P.A.  filaments  warm  up  at  least  3  minutes  before  the  high  voltage 
is  applied. 

Transmitter  peak  output  is  approximately  135  kW,  with  the  12.5  kV  plate  voltage 

3.5.  Power  Divider 

This  unit  divides  the  transmitter  power  output  into  two  parts,  A  and  B,  and 
selects  between  ordinary  (0)  and  extraordinary  (X)  modes  for  transmitting  with 
circular  polarization.  Mode  selection  is  determined  by  the  control  logic. 
During  0  transmission,  both  the  A  and  B  components  are  shifted  -90°  with 
respect  to  the  signal  from  the  P.A.  During  X  transmission,  the  shifts  are 
+90°  for  A  and  -90°  for  B.  See  Figure  3.3. 


8 


O  Polarization!  I  X  P 


transmitting  with  circular  polarization. 
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3.6.  Phase  Shifter 

The  phase  shifter  acts  only  on  the  A  component  from  the  power  divider.  For 
transmitting  in  either  the  0  or  X  modes  of  circular  polarization,  it  shifts 
the  A  component  by  another  -90°.  For  transmitting  with  linear  polarization, 
it  gives  no  shift  to  the  A  component.  See  Figure  3.3  for  details.  Mode 
selection  is  determined  by  the  control  logic. 

3.7.  Transmit  and  Receive  (T-R)  Switch 

The  T-R  switch  effectively  connects  the  antennas  in  the  array  to  the  transmit¬ 
ter  when  it  is  transmitting,  and  to  the  receiver  at  other  times.  All  signals 

thru  the  T-R  switch  are  shifted  in  phase  by  -90°.  As  noted,  the  T-R  must 
normally  pass  135  kW  of  power  during  transmission.  A  monitor  circuit  in  the 
T-R  circuit  provides  samples  of  the  A  and  B  components  for  monitoring  trans¬ 
mitter  performance.  Monitor  outputs  are  approximately  55  volts  peak  to  peak. 

3.8.  Matching  Network  Box 

There  are  two  matching-network  boxes,  one  for  each  of  the  two  T-R  switches. 

The  five  antenna  elements  aligned  in  the  east-west  direction  are  connected  to 
box  A,  and  the  five  north-south  elements  are  connected  to  box  B.  The  antennas 

with  impedance  50  ohms  each,  are  connected  in  parallel  to  the  outputs  of  the 

matching  network  boxes,  thus  providing  a  net  load  impedance  of  10  ohms  for 
each  of  the  two  boxes.  The  matching  network  matches,  in  both  cases,  the  50 
ohm  impedance  of  the  transmission  line  from  the  T-R  switch  to  the  10  ohm 
impedance  of  a  five-antenna  load. 

3.9.  Receiver  Polarization  Selector 

The  echo  signal  received  by  the  antennas  passes  thru  the  A  and  B  matching 
network  boxes  and  T-R  switches  to  the  receiver  polarization  selector.  The 
selector  has  two  input  ports,  A  and  B,  and  two  output  ports.  As  noted  above, 
the  phases  of  the  A  and  B  signals  are  both  shifted  by  -90°  in  the  T-R  switches 
If,  for  reception,  the  antenna  array  is  to  be  circularly  polarized,  the  re¬ 
ceiver  polarization  selector  adds  a  phase  shift  of  +90°  to  signal  A;  and  no 
shift  for  linear  polarization.  See  Figure  3.4. 
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Figure  3.4.  Phases  for  receiving  with  circular  polarization. 


The  output  ports  are  ports  2  and  4  of  a  4-port  hybrid  junction.  Signals  A  and 
B  are  applied  to  ports  3  and  1  of  the  hybrid,  respectively.  In  the  0  mode  of 
reception,  signals  at  ports  3  and  1  are  in  phase,  and  the  output  is  at  port  2. 
In  the  X  mode,  signals  at  ports  3  and  1  are  out  of  phase,  and  the  output  is  at 
port  4.  In  either  case,  the  non-used  output  port  (2  or  4)  is  terminated  with 
a  47  ohm  resistor.  Relays  K3  and  K4  select  the  appropriate  output  port. 
Control  of  the  receiving  polarization  mode,  circular  or  linear,  and  0  and  X 
selection  for  circular  polarization,  is  performed  by  the  control  logic. 

The  output  of  the  receiver  polarization  selector  passes  through  two  Hewlett- 
Packard  attenuators  in  series  that  can  be  set  over  the  ranges  0-12  dB  and 
0-120  dB,  respectively.  Appropriate  values  of  fixed  attenuation  are  selected 
during  normal  operation  and  during  system  calibration.  (See  Section  5  for 
details . ) 

3.10.  Receiver 


The  incoming  signal  from  the  two  fixed  attenuators  passes  through  a  series  of 
stepping  attenuators.  The  potential  sequence  of  attenuators  is  as  shown: 


Relays  determine  which  of  these  attenuators  is  in  the  series  at  any  moment, 
and  control  of  the  relays  is  performed  by  the  control  logic.  The  operator  can 
select  one  of  three  patterns  for  the  stepped  attenuation.  If  the  6  dB  pattern 
is  selected,  the  net  attenuation  of  the  series  is  stepped  through  values  of  0, 
6,  12,  and  18  dB.  If  the  9  dB  pattern  is  selected,  the  stepped  values  are  0, 
9,  18,  and  27  dB;  and  for  the  12  dB  pattern,  values  are  0,  12,  24,  and  36  dB. 
Output  of  the  stepping  attenuators  is  fed  to  the  rf  amplifier  of  the  receiver. 


The  rf  amplifier,  at  the  front  end  of  the  receiver,  has  a  low  or  high  gain 
mode,  switch  controlled.  The  gain  control  switch  sets  the  bias  on  the  rf 
amplifying  transistors.  The  gain  can  also  be  controlled  remotely,  via  a  cable 
connected  to  the  back  panel:  Terminating  the  cable  with  either  an  open  or  a 
short  sets  the  gain  remotely. 


22 


Receiver  i.f.  is  1.50  MHz,  and  the  local  oscillator  is  set  to  1.50  MHz  below 
the  rf  frequency  to  be  received.  The  i.f.  gain  can  be  controlled  by  a  trimpot 
located  on  the  front  panel.  Both  the  rf  and  i.f.  amplifiers  are  blanked  by 
the  PRF  control  signal  when  the  transmitter  is  on,  to  reduce  ringing  in  the 
receiver . 

Output  of  the  i.f.  amplifier  is  fed  to  a  square-law  detector  that  provides  an 
output  of  0  to  10  Vdc ,  in  proportion  to  the  power  of  the  received  signal. 

This  analog  output  passes  to  the  control  logic  where,  at  specified  sample 
times  (corresponding  to  sample  heights),  the  output  is  converted  to  a  6-bit 
digital  word. 

3.11.  Analog  To  Digital  Converter 

The  A/D  converter  is  located  on  card  8  of  the  control  logic  chassis.  It  is  a 
DATEL  converter  model  ADC-G8B1A,  with  8  bits  of  output.  However,  the  two 
least-significant  bits  are  not  used.  The  PDP-8  computer  has  a  12-bit  word 
size,  and  two  of  the  6-bit  words  from  the  converter  are  packed  into  one 
computer  word.  All  outputs  from  the  A/D  converter,  plus  corresponding  status 
information,  are  stored  (temporarily)  in  computer  memory. 

3.12.  Computer 

The  digital  computer  is  a  Digital  Equipment  Corporation  PDP-8/L  with  12-bit 
words,  parallel  input,  and  4,000  words  of  memory.  One  of  its  functions  is  to 
act  as  a  data  buffer  for  the  tape  recorder.  Another  function  is  to  average 
all  of  the  raw  data  in  a  run,  with  averages  for  each  sample  height  and  attenu¬ 
ation,  and  at  the  end  of  a  run,  print  the  results  on  the  teletype  for  on-site 
analysis.  These  averages  provide  a  field  record  of  D-region  condition  and 
equipment  operation.  (Further  explanations  of  the  computer  functions  are 
given  in  Section  5.) 

3.13.  Teletype 

The  ARS-33  teletype  prints  the  averages  calculated  by  the  computer,  after  a 
run  and  at  the  request  of  the  operator,  as  noted  in  Section  3.12.  It  can  also 
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print  the  contents  of  the  computer  memory  or  of  any  magnetic  tape  record,  for 
use  in  trouble  shooting  and  program  analysis,  or  for  verification  or  proper 
operation  of  the  recording  system. 

3.14.  Tape  Recorder  and  Controller 

The  magnetic  tape  recorder  is  a  Wang  Co.  model  10,  and  uses  a  lO1^  inch  reel  of 
computer  tape.  Recording  is  on  7  tracks  with  a  density  of  200  bpi.  The  Datum 
tape  controller  acts  as  interface  between  the  tape  recorder  and  the  computer. 
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SECTION  4 

A  n  tenna  Install at  ion 

4.1.  Antenna  Elements 


The  partial  reflection  sounder  has  been  operated  at  frequencies  2.2375  MHz  and 
2.6667  MHz.  The  2.6667  MHz  antenna  system  is  discussed  here,  and  similar 
concepts  apply  at  the  other  frequency. 


The  antenna  elements  are  made  of  copperweld  No.  10  wire.  Each  element  is  a 
center-fed  half-wave  dipole,  with  quarter-wave  arm  length  computed  from  the 
following  formula: 


3^iOj^ec_  *  0m  95  V  q  26.  7  ^6^ 

ic  jt  ^ 

(end-effect  -factor) 


#10  ju\  RE 
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Bottom  or  BAlu^) 
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Figure  4.1.  Dipole  Element. 


The  antenna  array  consists  of  five  crossed-d i pole  segments,  for  a  total  of  10 
half-wave  dipoles.  Figure  4.2  gives  the  dimensions  of  the  array. 


POLE'S  SO 
fve>OV/E  G>EOOMti 


(1)  Does  not  include  5’  for  pole  rope 

(2)  Includes  length  of  insulator 

(3)  Frequency  2.666666  MHz 

Figure  4.2.  Antenna  array,  plan  view. 


4.2.  Feed  Cables 


As  indicated  in  Figure  4.1,  a  feed  cable  and  a  balun  are  connected  to  the 
center  of  each  dipole.  Both  are  made  from  RG8  coaxial  cable.  Feed  cables 
connect  the  dipoles  to  the  matching  network  boxes,  and  are  one-half  wave,  with 
length  computed  as  follows: 
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In  practice,  the  feed  cables  are  cut  slightly  longer  than  the  computed  length, 

then  shortened  gradually  until  resonance  is  indicated  on  an  rf  bridge  (RX  meter, 

H.P.  model  250).  With  the  cable  unloaded,  typical  readings  at  resonance  on  the 

RX  meter  are  1.2  K  ohms  on  the  R  dial  and  ±5  pf  on  the  C  dial. 

P  P 


4.3.  Baiuns 


Baiun  length  is  computed  from  the  formula 


_A 
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As  with  the  feed  cables,  baluns  are  cut  slightly  longer  than  the  computed 
length,  then  shortened  to  resonance.  The  connection  of  the  balun  to  the 
feed  cable  is  shown  in  Figure  4.1.  The  feed  cable  and  the  balun  should 
be  taped  together  as  closely  as  possible  because  gaps  between  the  two  affect  the 
capacitive  reactance  of  the  antenna. 


*  The  factor  0.66  in  the  length  equations  is  applicable  to  RG8  cable. 
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4.4.  Antenna  Poles 


Thirteen  wooden  poles,  length  70  feet,  are  used  to  support  the  antennas. 

Figure  4.2  shows  the  distances  between  the  poles.  The  antennas  are  mounted  50 
feet  above  ground.  Once  the  antennas  have  been  calibrated,  their  height  must 
be  maintained  constant,  because  height  also  affects  the  capacitive  reactance. 

4.5.  Antenna  Calibration 

When  the  antennas  are  in  place,  the  connection  between  feed  cable  and  balun 
shields  is  made  approximately  one  quarter  wavelength  below  each  antenna.  The 
location  of  the  connection  is  then  varied  until  resonance  is  indicated  on  an 
RX  meter  attached  to  the  feed  cable.  At  resonance,  the  readings  should  be 
approximately  50  ohms  and  ±50  pf  . 

4.6.  Antenna  Matching 

After  all  antennas  are  calibrated,  the  two  matching  network  boxes  must  be 
calibrated  as  shown  in  Figure  4.3  (for  box  A). 


RX 

MAvTCH 

NETWORK 

METER 

R&8 

Box 

A 

~  -feed  cables  to 

east -west  dipoles 


Figure  4.3.  Antenna  matching  calibration  setup. 

The  RX  meter  is  set  at  50  ohms  R  and  0  pf  C  ,  and  connected  to  the  matching 

P  F  P 

network  box  with  A  calibrated  at  50  ohms  half-wavelength  cable.  The  box  is 
then  adjusted  to  resonance  using  the  variable  capacitor  and  inductor  within 
the  box. 
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SECTION  5 


System  Operation 

5.1.  Transmitter  Turn-on 


To  activate  the  transmitter  and  associated  equipment,  turn  on  the  following 
in  sequence: 

a.  Main  power  (SI) 

b.  Exciter  high  voltage  (S2) 

c.  P.  A.  filaments  (S3);  adjust  to  7.4  Vac 

d.  Pulse  generator 

e.  Receiver 

f.  P.  A.  high  voltage  (S4) ;  adjust  to  12.5  kVdc 


5.2.  Logic  Turn-on 

To  activate  the  logic  unit  and  associated  equipment,  turn  on  the  following 
in  sequence: 


a.  Relay  power  supply 

b.  Logic  unit  power  supply 

c.  Tape  controller 

d.  Tape  recorder 

e.  Computer 

f.  Teletype 


5.3.  Program  Loading 

To  load  the  operating  program  into  the  computer,  do  the  following: 

a.  Thread  the  program  tape  onto  the  tape  drive  and  press  the 

"load"  switch.  The  tape  should  advance  to  the  start  marker, 
and  the  "on  line"  status  light  should  come  on. 
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b.  With  the  computer  in  single-step  mode,  load  address  7600,  which 
is  the  entry  address  for  the  bootstrap  loader.  The  contents  at 
address  7600  will  be  displayed.  Step  thru  and  examine  the 
remainder  of  the  bootstrap,  shown  in  Table  5.1.  To  see  the 
contents  of  the  second  address,  7601,  press  the  ’’examine’1 
switch  twice.  Each  subsequent  address  can  then  be  seen  via 
one  operation  of  the  "examine”  switch.  If  any  parts  of  the 
bootstrap  program  are  not  correct,  load  the  correct  instruc¬ 
tions  manually.  This  is  done  by  re-loading  the  address  in 
question,  loading  the  correct  contents  on  the  registers,  and 
then  pressing  the  "deposit"  switch. 

c.  Go  back  to  address  7600  and  deactivate  the  single-step  mode. 

d.  Clear  the  registers  (0000)  and  press  the  "start"  switch. 
Execution  will  halt  at  address  7603. 

e.  Load  instruction  0420  on  the  registers  and  press  the  "continue" 
switch.  The  accumulator  will  then  indicate  0001.  If  this  is 
not  the  case,  examine  the  bootstrap  procedure  and  try  the 
previous  steps  again.  The  accumulator  must  read  0001  if  the 
operating  program  is  to  function  properly. 

5.4.  Data  Sampling  Procedure 

a.  Load  a  clean  data  tape  and  make  sure  the  handler  is  "on  line."  An 
end-of-file  (EOF)  mark  must  be  entered  at  the  start  of  the  tape  and 
after  each  data  run.  To  write  the  EOF  mark,  enter  instruction  4000 
at  address  0200,  then  press  the  "continue"  switch.  The  tape  should 
advance  as  the  EOF  is  written. 

b.  Enter  instruction  0200  (initialize)  at  address  0200,  then  press  the 
"start"  switch.  The  teletype  will  print  "CN="  and  the  operator  then 
types  "1.0  return. "  This  indicates  the  ratio  between  0  and  X  ampli¬ 
tudes  . 
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TABLE  5.1.  COMPUTER  BOOTSTRAP  PROCEDURE 


Address 

Contents 

7600 

6711 

7601 

5200 

7602 

6706 

7603 

7402 

7604 

7300 

7605 

7604 

7606 

6716 

7607 

6722 

7610 

7300 

7611 

6713 

7612 

5211 

7613 

6726 

7614 

3623 

7615 

2223 

7616 

2222 

7617 

5211 

7620 

6723 

7621 

5200 

*  7622 

0177 

*  7623 

0000 

*  After  the  operating  program  has  been  loaded,  the  contents  of 
address  7622  and  7623  will  be  0000  and  7601,  respectively. 
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c.  The  fixed  attenuators  at  the  receiver  input  should  now  be  adjusted 
for  proper  attenuation.  This  requires  that  the  transmitter  be  in 
operation.  Press  the  "start"  button  on  the  control  logic  unit.  The 
transmitter  will  start  to  operate,  but  the  received  signals  will  not 
be  saved  either  in  the  computer  or  on  the  data  type.  The  proper 
amount  of  fixed  attenuation  is  determined  by  visual  inspection  of 
the  receiver  oscilloscope  display.  Normal  fixed  attenuation  for 
daytime  operation  is  the  range  6  to  10  dB.  To  stop  transmitting,  it 
is  necessary  to  reset  the  control  logic.  This  is  done  by  pushing 
the  button  switch  on  card  3  of  the  control  unit,  located  behind  the 
front  panel. 

d.  Enter  the  appropriate  date,  time,  etc.,  via  the  thumbwheels  on  the 
front  panel  of  the  control  unit.  This  information  will  be  recorded 
on  the  tape,  along  with  the  data. 

e.  Clear  (0000)  the  computer  registers  and  press  the  computer  "continue" 
switch.  The  computer  is  now  ready  to  receive  the  data. 

f.  Press  the  "start"  button  on  the  control  logic  unit.  This  starts  the 
run,  with  the  transmitter  operating,  and  the  received  signals  digi¬ 
tized.  The  computer  will  store  the  digitized  data,  and,  periodically, 
dump  the  data  onto  the  tape.  The  time  interval  for  a  normal  data 
run  is  10  to  15  minutes. 

g.  To  terminate  the  run,  press  the  "stop"  button  on  the  control  unit. 
Then  write  an  EOF  mark  on  the  data  tape.  (See  step  5.4a.) 

h.  (Optional)  To  print  out  the  data  averages  on  the  teletype,  for  the 
run  just  concluded,  load  instruction  2000  ("computations")  at  address 
0200  and  press  the  "continue"  switch.  To  stop  computations,  press 
the  "stop"  switch. 
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5.5.  Receiver  Calibration 


The  purpose  of  the  receiver  calibration  is  to  determine  the  linearity  and 
sensitivity  of  the  receiver,  and  the  proper  functioning  of  the  stepped 
attenuators.  Figure  5.1  shows  the  calibration  setup. 


H  P.  -fixed 
attenuators 


Figure  5.1.  Receiver  Calibration  Setup. 

Normal  saturation  threshold  for  the  receiver  is  approximately  -97  dBm  at 
zero  attenuation.  Normal  calibration  levels  are  as  follows:  -90,  -96, 
-98,  -100,  -102,  and  -108  dBm.  These  levels  are  selected  using  the  fixed 
attenuators.  At  each  level,  record  a  "run"  of  (calibration)  data  on  the 
tape,  of  several  minutes  duration.  Note  that,  in  addition  to  the  fixed 
attenuation,  the  stepped  attenuators  will  be  stepping  through  their  four 
levels  during  the  run.  Use  a  step  unit  of  6dB.  Calibration  data  is 
written  on  the  data  tape  in  exactly  the  same  format  as  regular  data.  A 
printout  of  averages  on  the  teletype  shows  the  operation  of  the  stepping 
attenuators  and  the  receiver  linearity. 
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5.6.  Turn-off  Procedure 

To  shut  down  the  station  turn  off  the  following  in  sequence: 

a.  Teletype 

b.  Computer 

c.  Tape  drive 

d.  Tape  controller 

e.  Logic  unit  power  supply 

f.  Relay  power  supply 

g.  Receiver 

h.  Pulse  generator 

i.  Reduce  P.  A.  high  voltage  to  zero,  then  turn  off  (S4) 

j.  Exciter  high  voltage  (S2) 

k.  Wait  5  minutes  for  the  P.  A.  tubes  to  cool,  then  turn  off 
the  P.A.  filaments  (S3) 

l.  Main  power  (SI) 
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SECTION  6 


Data  Processing 

The  first  step  in  the  data  processing  is  to  time-average  the  data  recorded  on 
magnetic  tape  during  a  run.  This  averaging  is  done  by  a  computer  program  that 
runs  at  White  Sands  Missile  Range  on  the  UNIVAC  1108  system.  A  listing  of  the 
program,  which  is  written  in  FORTRAN,  is  given  in  Appendix  II.  Our  intention 
here  is  to  summarize  what  this  program  does. 

Recall  that  each  pulse  from  the  transmitter  results  in  an  echo  that  is  digi¬ 
tized  via  30  samples,  and  that  8  consecutive  pulses  form  the  basic  operating 
pattern.  Thus  in  averaging  the  data,  240  separate  averages  must  be  computed. 
We  actually  perform  480  averages  because  of  the  fact  that  the  raw  data  is 
passed  through  two  different  screening  procedures.  (Which  of  the  two  results 
will  be  used  is  decided  after  the  fact.)  A  slight  complication  arises  in  that 
each  record  on  the  data  tape  contains  two  patterns,  i.e.,  echoes  from  16 
consecutive  pulses,  and  thus  yields  two  numbers  for  each  running  average. 

Now  each  data-number  is  an  integer  between  0  and  63,  and  its  value  is  a  mea¬ 
sure  of  received  signal  strength.  "Strength"  is  the  equivalent  of  signal 
intensity;  but  the  information  required  is  signal  amplitude .  Each  data-number 
is  converted,  prior  to  averaging,  into  a  corresponding  measure  of  signal 
amplitude  via  a  table-lookup  procedure  that  is  based  on  receiver  calibration. 

Appendix  IV  contains  details  regarding  receiver  calibration.  We  shall  simply 
state  here  that  the  end-product  of  receiver  calibration  is  a  table  or  64 
numbers  scaled  so  that  the  first  is  0  and  the  last  is  63.  (Scaling  to  this 
range  is  done  only  for  convenience.  Since  the  end  result  is  the  amplitude 
ratio  Ax/Aq,  any  constant  used  to  scale  the  amplitudes  cancels  out  in  the 
ratio.)  This  table  is  named  RX  in  the  program.  If,  for  example,  the  consi¬ 
dered  data-number  is  15,  then  the  corresponding  measure  of  signal  amplitude  is 
the  16th  value  in  the  RX  array,  RX(16). 
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The  "heart"  of  the  averaging  program  is  the  accumulation  loop,  reproduced  in 
part  in  Figure  6.1.  Therein,  the  accumulations  are  carried  out  for  each  of 
the  separate  averaging  calculations.  Comments  follow  for  each  of  the  state¬ 
ments  identified  by  number  in  Figure  6.1: 


1  DO  210  J=1 ,8 

2  K=J 

3  195  Kl=l 

4  K2=l 

5  IF  (IDATA(K,IJKHT) .GT.MAX1)  K1=0 

6  N1 (J)=N1 (J)+K1 

7  IF  (IDATA(K,IJKHT) .GT.MAX2)  K2=0 

8  N2(J)=N2(J)+K2 

9  DO  205  H=1 ,MM 

10  IF  (NSEG.EQ. 1)  GOTO  197 

11  IF  (IDATA(K,H) .GT.MAXX)  N3(J,H)=N3(J,H)+1 

12  I=IDATA(K,H)+1 

13  A1(J,H)=A1(J,H)+K1*RX(1) 

14  A2(J,H)=A2(J,H)+K2*RX(1) 

15  GOTO  200 

16  197  A1 (J,H)=IDATA(J,H) 

17  A2(J,H)=IDATA(J+8,H) 

18  200  CONTINUE 

19  205  CONTINUE 

20  K=K+8 

21  IF  ((K.LE. 16) .AND. (NSEG.NE. 1))  GOTO  195 

22  210  CONTINUE 


Figure  6.1.  Accumulation  logic  for  data  averaging. 

These  statements  are  from  the  computer 
program  in  Appendix  II. 
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1.  J  (from  1  to  8)  indexes  the  pulses  of  the  basic  8-pulse  pattern. 

2.  K  is  initially  the  same  as  J. 

3,4.  Multipliers  K1  and  K2  are  initially  set  at  1. 

5.  IDATA  (size  16  by  30)  contains  the  480  data-values  from  one  record 

of  the  data  tape;  16  echoes,  30  samples  per  echo.  Index  IJKHT  is  a 
constant  set  by  card  input.  If,  for  example,  IJKHT=4,  then  K1  is 
changed  to  zero  if  IDATA(K,4)  is  greater  than  MAXI,  which  is  another 
constant  set  by  card  input.  The  idea  here  is  this:  assume  that  the 
return  from  the  4th  (if  IJKHT=4)  sampled  altitude  is  sure  to  contain 
noise  only,  no  partial  reflections  being  expected  from  this  particu¬ 
lar  altitude.  Then  if  the  noise  (the  data  value)  here  exceeds  the 
threshold  value  MAXI,  all  30  values  from  this  echo  are  not  going  to 
be  included  in  the  averages  (thanks  to  setting  K1=0). 

6.  N1(J)  counts  the  number  of  times  that  K1  remains  at  1,  which  corre¬ 
sponds  to  the  number  of  values  averaged. 

7.  Screening  similar  to  that  at  line  5  is  done  here,  but  at  a  different 
threshold,  MAX2 . 

8.  Counting  similar  to  that  at  line  6. 

9.  H  (from  1  to  MM=30)  indexes  the  samples  in  each  echo. 

10.  NSEG  is  a  contant  set  by  card  input.  It  fixes  the  number  of  records 
to  be  accumulated  prior  to  averaging.  If  all  records  in  a  run  are 
to  be  included,  any  value  for  NSEG  larger  than  the  number  of  records 
will  do.  Say  a  run  contains  800  records.  If  NSEG=100,  then  the  run 
is  effectively  broken  into  8  parts,  and  averaging  will  be  done  over 
each  of  these  parts,  in  turn.  NSEG=1  is  a  special  case  that  results 
in  the  listing  of  all  data  values  recorded  on  the  tape  for  this  run 
(480  values  per  record,  and  one  record  per  page  of  the  printout). 

11.  If  data  value  IDATA(K,H)  exceeds  a  value  MAXX,  counter  N3(J,H)  is 
incremented.  We  use  MAXX=62  so  that  N3  counts  the  number  of  samples 
for  which  the  digitizer  was  probably  saturated.  Any  averages  for 
which  the  N3  count  is  not  small  are  suspect. 

12.  The  data  value  is  incremented  by  one  and  the  result,  I,  will  serve 
as  an  index  for  the  receiver-calibration  array  RX  (size  64). 

13.  Array  A1  (size  8  by  30)  contains  the  240  accumulations  corresponding 
to  screening  via  MAX  1.  RX(I)  is  the  data  value  that  will  be  ac¬ 
cumulated  if  K1  is  not  zero. 
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Array  A2  contains  the  240  accumulations  corresponding  to  screening 
via  MAX2. 

16,17.  In  the  special  case  NSEG=1,  the  first  240  values  of  the  record  are 
copied  into  A1  (without  conversion  via  RX) ,  and  the  second  240 
values  are  copied  into  A2.  The  A1  and  A2  arrays  will  be  printed 
before  the  next  record  is  processed. 

19.  End  of  inner  loop  on  sample  (altitude)  index  H. 

20.  Increase  K  by  8. 

21.  If  NSEG>1  and  if  K^16,  go  back  and  accumulate  values  from  the  second 
half  of  the  record  (echoes  9  through  16). 

22.  End  of  outer  loop  on  echo  (pulse)  index  J. 

An  example  of  the  output  obtained  from  the  averaging  program  is  shown  in 

Figure  6.2.  A  page  of  output  similar  to  this  is  generated  for  every  run,  or 

for  every  run  segment,  that  is  averaged. 

The  following  comments  relate  to  Figure  6.2,  and  describe  the  material  found 

on  the  lines  marked,  at  the  left  of  the  figure,  with  letters  A  through  F: 

A.  •  Tape  file  number. 

•  Starting  height;  i.e.,  the  lower  limit  of  the  altitude  interval 
sampled . 

•  Segment  number  (applicable  to  averages  over  parts  of  the  run). 

•  Records  per  segment  (NSEG) .  Since  999  is  greater  than  the 
number  of  records  here,  the  entire  run  is  processed  as  one 
segment. 

•  Reference  height  index  (IJKHR)  used  to  select  noise  values  from 
each  echo  to  test  against  the  screening  thresholds  MAXI  and 
MAX2 . 

B.  •  Time,  date,  and  comments  —  reference  information  entered  via 

cards . 

•  Solar  zenith  angle,  computed  using  the  preceeding  time  and 
date.  See  Appendix  III. 

C.  Reference  information  contained  on  the  data  tape,  and  obtained 
from  thumbwheels  set  by  the  operator  prior  to  the  run. 
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Number  of  records  processed  (in  this  case,  all  records  in  the 
file) . 

Selected  screening  threshold  values;  10  for  MAX  1  and  5  for 
MAX2 . 

Parity  check  .  .  .  The  number  of  records  read  from  the  tape 
without  parity  error.  Since  this  number  here  equals  the  number 
of  records  used,  there  were  no  parity  errors. 

Reference  information  from  the  tape  identifying  the  particular 
mode  of  system  operation  (normal  operating  mode,  in  this  case). 
Average  amplitudes,  arranged  in  16  columns,  with  rows  holding 
results  from  55  km  to  113  km  (as  per  the  left-most  column). 
Columns  1  through  8  (see  numbers  below  the  columns)  contain  the 
240  averages  for  MAXI  screening,  while  columns  9  through  16 
contain  the  averages  for  MAX2  screening.  Column-sets  1-8  and 
9-16  each  follow  the  basic  8-pulse  pattern:  amplitudes  for 
ordinary  mode  with  0,  1,  2,  and  3  units  of  attenuation;  fol¬ 
lowed  by  amplitudes  for  extraordinary  mode  with  0,  1,  2,  and  3 
units  of  attenuation.  (The  attenuation  unit  here  is  6  dB.) 

There  is,  in  the  receiver,  a  delay  corresponding  to  an  altitude  increment  of 
5  km.  This  delay  is  not  taken  into  account  in  the  timing  circuits  that  set 
(via  operator  selection)  the  starting  height.  The  delay  ijs  taken  into  account 
in  the  data  processing  that  follows.  Thus  the  data  in  Figure  6.2  spans  a  true 
altitude  interval  of  50  km  to  108  km,  rather  than  the  indicated  interval  55  km 
to  113  km.  The  operator  takes  note  of  this  delay.  For  example,  if  he  wants 
to  start  sampling  at  50  km,  he  will  select  a  value  indicated  as  55  km  on  his 
control  unit.  We  have  found  it  less  confusing  to  use  the  indicated  value 
within  the  averaging  program. 

Counts  accumulated  by  N1(J),  N2(J),  and  N3(J,H)  (see  Figure  6.1)  are  printed 
on  a  separate  page  of  computer  output.  For  example,  in  the  run  of  Figure  6.2, 
N1  showed  that  1035,  1466,  1726,  2744,  1099,  1436,  1719,  and  1744  values  were 
included  in  the  averages  for  columns  1-8,  respectively.  AU  data  was  thus 
used  in  the  averages  for  columns  4  and  8,  as  indicated  by  N 1  =  1 744=2 '872 ,  872 
records  having  been  considered.  The  N3  counter  showed  that  only  a  relatively 
few  saturated  values  were  included  in  the  averages  of  columns  3,  4,  7,  and  8, 


D.  • 

E. 

F. 
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and  only  at  some  upper  altitudes.  Columns  1,  2,  5,  and  6,  however,  contain  an 
excessive  number  of  saturated  values  at  many  altitudes,  and  should  not  be 
used . 

For  any  run,  the  amplitudes  AQ(h)  and  A^(h)  obtained  from  the  averaging  pro¬ 
gram  (Appendix  II)  are  used  to  obtain  the  electron  density  profile  N(h)  for 
the  run.  The  procedure  here  requires  the  amplitude  ratio  A^/A^,  to  be  used  as 
in  Equation  (1.5). 

The  first  step  is  to  compute  by  hand  an  A^/A^  vs.  altitude  profile  using 
average  amplitude  values  such  as  found  in  Figure  6.2.  Because  of  the  fact 
that  data  were  recorded  at  four  different  values  of  receiver  attenuation,  and 
that  raw  data  were  screened  using  both  MAXI  and  MAX2  thresholds,  the  analyst 
must  decide  which  of  several  possible  values  are  to  be  used  to  compute  A^/A^ 
at  each  altitude.  Judgements  here  must  be  based  on  experience  with  similar 
data.  Because  of  the  subjective  nature  of  the  decision  process,  we  have  not 
attempted  to  automate  this  step. 

The  following  comments  are  based  on  the  example  run  of  Figure  6.2,  and  are 
intended  to  provide  some  insight: 

•  "Amplitudes"  below  67  km  (indicated  altitude)  are  simply  noise,  and 
contain  no  useful  information.  In  fact,  we  use  no  data  below  75  km. 
As  noted,  an  indicated  altitude  of  75  km  on  the  listing  represents  a 
true  altitude  of  70  km,  etc. 

•  A^  values  at  105  km  (listed  altitude)  and  above  are  simply  noise  and 
cannot  be  used.  In  fact,  we  use  no  data  above  89  km  (84  km  true). 

•  There  is  little  difference  between  the  averages  in  columns  1-8  and 
9-16,  showing  that  MAXI  and  MAX2  thresholds  are  nearly  equivalent. 

But  since  more  values  are  included,  generally,  under  the  MAXI  screen¬ 
ing,  we  will  use  only  columns  1-8. 

•  The  "best"  data  seems  to  be  in  columns  3  and  7.  Amplitudes  in 
columns  1,  2,  5,  and  6  tend  to  be  too  large,  indicating  receiver 
saturation,  while  those  in  columns  4  and  8  tend  to  be  too  small. 
(Receiver  attenuation  is  stepped  so  that  choices  such  as  this  a  re 
possible ! ) 
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The  A^/A^  calculation  then  proceeds  as  shown  in  Table  6.1: 


TABLE  6.1.  AMPLITUDE  CALCULATIONS  FOR  DATA  IN  FIGURE 


true 

altitude,  km  A  A  A  /A 

*  0  X  X  o 


70 

6 

732 

10 

371 

1 

54 

72 

14 

724 

19 

832 

1 

35 

74 

20 

925 

24 

783 

1 

18 

76 

25 

767 

28 

617 

1 

11 

78 

28, 

.554 

29. 

.391 

1, 

.03 

80 

29. 

.963 

20. 

.425 

0, 

.68 

82 

32. 

.877 

17, 

.700 

0, 

.42 

84 

31. 

.699 

13 

.152 

0 

.41 

6.2. 
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The  data  used  in  obtaining  the  averages  in  columns  3  and  7,  Figure  6.2,  were 
all  obtained  at  the  same  receiver  attenuation,  and  so  the  averages  can  be  com¬ 
bined  directly  to  obtain  A^/A^.  If  one  wishes  to  use  values  from  non-paired 
columns,  compensation  must  be  made.  For  example,  say  we  choose  to  use  A^  = 
12.401  from  column  8  at  74  km  (true)  rather  than  24.783  from  column  7.  Then 
A^/A^  =  12.401/20.925  =  0.593.  For  this  run,  the  level  of  stepped  attenuation 
was  set  at  6dB,  so  the  amplitudes  in  column  8  are  low,  relative  to  columns  3 
and  7,  by  6dB,  or  a  factor  of  1.995.  Thus,  the  preceding  value  of  0.593  must 

be  increased  by  this  factor:  A  /A  =  0.593  •  1.995  =  1.18. 

J  x  o 

Accordingly,  one  might  expect  all  values  in  column  7,  say,  to  be  essentially 
twice  the  corresponding  values  in  column  8.  This  relationship  does  hold, 
approximately,  for  mid-value  amplitudes,  but  is  clearly  not  the  case  for  low 
amplitudes  (due  to  noise  effects)  and  high  amplitudes  (due  to  saturation 
effects) . 

The  A  /A  profile  in  Table  6.1  above  is  typical  in  that,  with  increasing  alti¬ 
tude,  it  drops  quite  smoothly  from  values  above  1  to  values  below  1.  In  some 
runs,  the  ratio  falls  to  values  of  0.1  or  less  at  high  altitude,  rather  than 
the  limit  of  about  0.4  in  this  run.  Also,  some  profiles  show  a  maximum  near 
the  low-altitude  end.  In  this  run,  if  a  maximum  does  exist,  it  falls  short  of 
the  range  of  useful  data.  It  should  be  noted  that  our  decision  to  limit  the 
data  interval  to  between  70  km  and  84  km  was  based  primarily  on  the  A^/A^ 
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profile,  as  results  outside  this  interval  seem  unrealistic  (i.e.,  non  physical). 
Our  computer  program  that  changes  the  8  values  of  A^/A^  in  Table  6.1  (1.54  to 
0.41)  into  an  electron  density  profile  is  shown  in  Figure  6.3.  The  program 
name  is  RLTIM3  and  the  language  is  API.  As  per  Equation  (1.5),  the  program 
must  access  tables  R(h)  and  G(h) .  Examples  of  such  tables  are  shown  in  Table 
6.2,  with  Version  C  appropriate  to  the  run  we  are  considering.  The  particular 
tables  to  be  used  by  RLTIM3  are  to  be  named  RTEST  and  GTEST. 

The  call  of  RLTIM3  includes  arguments  A  and  AXAO.  AXAO  is  the  list  of  values 
of  A^/A^,  namely  the  8  values  from  Table  6.1  for  the  example  at  hand.  Argument 
A  consists  of  3  integers: 

•  A(l)  is  the  altitude  in  kilometers  of  the  first  entry  in  AXAO  (70 
here) ; 

•  A(2)  is  the  altitude  of  the  first  value  of  AXAO  to  be  utilized 
(again,  70  here) ; 

•  A(3)  is  the  degree  of  the  polynomial  to  be  fit  to  the  values  of 

y (h)  =  In  R/AXAO 

(y  values  are  computed  at  line  8  of  the  program).  In  our  example 
run,  a  value  A(3)  =  4  was  found  to  provide  good  results.  Generally, 
the  value  should  be  about  equal  to  half  the  number  of  data  points  (8 
here) . 

At  line  11  of  the  program,  a  procedure  GO*  is  called  to  perform  the  least- 
squares  eval'jition  of  the  polynomial  coefficients  a^.  These  coefficients 
provide  "best"  values  (y)  of  y  at  each  altitude  (x)  over  the  range  considered: 

2  3 

y  =  a^  +  a^x  +  a^x  +  a^x  . 

Although  the  altitude  interval  between  data  samples  is  2  km,  we  compute  results 
at  1  km  intervals  using  the  preceding  equation. 


*G0,  in  turn,  calls  procedure  ITER  which  calls  procedure  FUN.  Listings  of 
GO,  ITER,  and  FUN  are  included  in  Figure  6.3. 
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VRLTIM 3[[]]7 

V  A  RLTIM3  AXAO 

[I]  SEQ*-  012345 
C  2  ]  N+pAXAO 

[3]  NN*-p((  {AZ21-AZ1'\)*2)*AXA0) 

[4]  tf0«-/l[2] 

[  5  ]  N1+-A11] 

[6]  N3*-A [  3  ] 

[7]  N2<-NU-2*(N-1) 

[8]  Y-—RTESTl(ltO-SO)  +  2xO.\(NN-l)lH  (A  [  2  ] -A  C  1] )  *  2  )  *AXAO 

[9]  jr«-2xo,t2x(flW-l)*2 

CIO]  A*-N3p0 

[II]  GO 

[12]  A-*-6pA  ,0, 0,0,0 

[13]  X*-0  ,  i  2*NN-1 

[14]  EDEN4-At2]  +  (2*AL3l*X)  +  (  3xAtHlxXxX)  +  (t*xA[5lxX*XxX)  +  5xAlCy'\xXxXxXxX 

[15]  EDEN*-EDENiGTESTl (W0-50 )+X] 

[16]  IEDEN+0.5 

V 

VG0[[]]7 

7  (7(7 

[1]  a  LSQ  CURVE  FIT  PROCEDURE 

[2]  COUNT*- 1 

[3]  FJ1  s  ( pi4  )  irmpF) 

[4]  C0UNT+-C0UNT+ 1 

[5]  +(C0UNT<3)/FIl 

[  6  ]  £7?P 

[7]  0. 01x[0. 5+100xFF 

V 

Vim?[[]]7 
7  K  ITER  N  ;P 


[1] 

a  LSQ  CURVE  FIT  PROCEDURE 

[2] 

FF*-  OpO 

[3] 

J+ERR+0 

[4] 

P*-(K,K)p0 

VFUN cniv 

[5] 

l)p0 

7  A  FUN  X;N 

[6] 

iri:W  +  l 

[1] 

n+pA 

[7] 

A  FUN  XlJ] 

[2] 

V+6p0 

[8] 

F+G-TlJl 

[3] 

4«-6p;4  f  0 , 0 , 0 , 0 

[9] 

FF*-FF  ,  F 

[4] 

V<rX*$EQ 

[10] 

ERR*-ERR+F*F 

[5] 

G+-  +  /  A  x  V 

[11] 

DH*-{  \,K)pV 

[6] 

V+NpV 

[12] 

DV*-{K,\)pV 

[7] 

A+NpA 

[13] 

P+P+DV+ .*PR 

T7 

[14] 

M*-M+F*DV 

[15] 

■+(J<N)/IT1 

[16] 

P«-SP 

[17] 

AA+-P+  .  *M 

[18] 

AA+KpAA 

[19] 

A+A-AA 

[20] 

ERR*-(.ERRi(pX)-(pA)  )*0. 5 

7 


Figure  6.  3.  Procedures  used  in  standard  data  processing. 
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TABLE  6.2.  TABULATED  FUNCTIONS  R(h)  AND  G(h)  USED  IN  THE  STANDARD  DATA 
PROCESSING.  VERSIONS  A,  B,  AND  C  WERE  COMPUTED  USING  THE 
CONDITIONS  DEFINED  IN  TABLE  6.3 


Altitude  R(h)  G(h) 


h,  km  A  B 


51 

1.1458 

1.1435 

52 

1.1679 

1 . 1645 

53 

1.1880 

1.1835 

54 

1.2135 

1.2074 

55 

i  .2394 

1.2315 

56 

1.2733 

1.2628 

57 

1.3085 

1.2950 

58 

1.3472 

1.3300 

59 

1.3904 

1.3685 

60 

1.4378 

1.4103 

61 

1.4866 

1.4526 

62 

1.5453 

1.5027 

63 

1.6107 

1.5575 

64 

1.6816 

1.6159 

65 

1.7594 

1.6786 

66 

1 . 8442 

1.7456 

67 

1.9393 

1.8191 

68 

2.0457 

1.8990 

69 

2.1574 

1.9805 

70 

2.2731 

2.0621 

71 

2.3908 

2.1421 

72 

2.5103 

2.2202 

73 

2.6323 

2.2963 

74 

2.7546 

2.3692 

75 

2.8757 

2.4376 

76 

2.9949 

2.5014 

77 

3.1077 

2.5586 

78 

3.2128 

2.6091 

79 

3.3044 

2.6507 

80 

3.3850 

2.6855 

81 

3.4513 

2.7128 

82 

3.5055 

2.7340 

83 

3.5477 

2.7499 

84 

3.5801 

2.7615 

85 

3.6048 

2.7701 

86 

3.6229 

2.7762 

87 

3.6361 

2.7804 

88 

3.6456 

2.7834 

89 

3.6522 

2.7855 

90 

3.6569 

2.7869 

91 

3.6603 

2.7879 

92 

3.6626 

2.7886 

93 

3.6643 

2.7891 

94 

3.6654 

2.7895 

95 

3 . 6663 

2.7897 

96 

3.6668 

2.7899 

97 

3.6672 

2.7900 

98 

3.6674 

2.7901 

99 

3.6676 

2.7901 

100 

3.667 

2.7902 

C  A 

1.1933  2 . 2568E-5 
1.2223  2 . 9385E-5 
1.2488  3 . 6247E-5 
1.2823  4.5710E-5 
1.3162  5.6109E-5 
1.3607  7 . 0852E-5 
1.4069  8 . 7286E-5 
1.4575  1.0647E-4 
1.5137  1.2894E-4 
1.5754  1.5471E-4 
1.6386  1 . 8198E-4 
1.7145  2 . 1548E-4 
1.7987  2.5310E-4 
1.8897  2.9368E-4 
1.9890  3 . 3720E-4 
2.0967  3.8287E-4 
2.2165  4.3108E-4 
2.3491  4 . 8052E-4 
2.4865  5 . 2662E-4 
2.6270  5 . 6756E-4 
2.7673  6.0155E-4 
2.9071  6 . 2773E-4 
3.0468  6. 4530E-4 
3.1835  6 . 5303E-4 
3.3153  6.5012E-4 
3.4414  6 . 3585E-4 
3.5573  6 . 1042E-4 
3.6620  5 . 7437E-4 
3.7505  5 . 3106E-4 
3.8262  4.8139E-4 
3.8867  4 . 2984E-4 
3.9350  3. 7799E-4 
3.9716  3.2906E-4 
3.9991  2 . 8398E-4 
4.0196  2 . 4295E-4 
4.0343  2.071  IE-4 
4.0448  1 . 7599E-4 
4.0523  1 . 4908E-4 
4.0574  1 . 2638E-4 
4.0610  1.0715E-4 
4.0636  9 . 0626E-5 
4.0654  7 . 6561E-5 
4.0666  6 . 4433E-5 
4.0675  5.3816E-5 
4.0681  4 . 4792E-5 
4.0685  3 . 7234E-5 
4.0688  3 . 0923E-5 
4.0690  2 . 6192E-5 
4.0691  2 . 1434E-5 
4.0692  1 . 7863E-5 


B  C 

2 . 1942E-5  2.8810E-5 
2 . 8295E-5  3.7186E-5 
3 . 4608E-5  4 . 5523E-5 
4 . 3196E-5  5 . 6886E-5 
5 . 2492E-5  6.9213E-5 
6 . 5440E-5  8 . 6435E-5 
7 . 9581E-5  1 . 0532E-4 
9.5726E-5  1.2698E-4 
1 . 1418E-4  1.5189E-4 
1 . 3478E-4  1 . 7990E-4 
1 . 5594E-4  2.0892E-4 
1 . 8107E-4  2 . 4376E-4 
2 . 0820E-4  2.8187E-4 
2.3614E-4  3 . 2180E-4 
2 . 6458E-4  3.6327E-4 
2 . 9264E-4  4.0552E-4 
3.2012E-4  4 . 4763E-4 
3 . 4566E-4  4.8886E-4 
3 . 6654E-4  5.2478E-4 
3.8186E-4  5 . 5389E-4 
3.9102E-4  5.7499E-4 
3.9393E-4  5.8765E-4 
3. 9049E-4  5.9142E-4 
3. 8071E-4  5 . 8580E-4 
3.6490E-4  5.7071E-4 
3 . 4334E-4  5.4606E-4 
3.1727E-4  5 . 1302E-4 
2 . 8766E-4  4.7265E-4 
2.5714E-4  4 . 2865E-4 
2 . 2600E-4  3.8170E-4 
1 . 9655E-4  3 . 3563E-4 
1.6904E-4  2.9131E-4 
1 . 4456E-4  2.5094E-4 
1 . 2302E-4  2.1475E-4 
1.0410E-4  1 . 8252E-4 
8 . 8022E-5  1 . 5482E-4 
7 . 4341E-5  1 . 3107E-4 
6 . 2687E-5  1.1072E-4 
5 . 2969E-5  9.3672E-5 
4.4801E-5  7 . 9298E-5 
3.7827E-5  6.6999F.-5 
3.1916E-5  5.6557E-5 
2 . 6835E-5  4.7570E-5 
2 . 2398E-5  3.9715E-5 
1 . 8633E-5  3. 3045E-5 
1 . 5484E-5  2.7464E-5 
1 . 2857E-5  2 . 2806F.-5 
1 . 0888E-5  1.9315E-5 
8 . 9093E-5  1 . 5805E-5 
7 . 4247E-5  1.3172E-5 
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TABLE  6.3.  SPECIFICATIONS  USED  IN  COMPUTING  THE  FUNCTIONS  OF  TABLE  6.2 


All  functions  were  computed  using  the  Flood  approximation,  Equations  1.3A 
and  1.4A,  and  the  collision  frequency  profile  from  Table  7.1. 


VERSION 

A 

B 

C 

For  station  located  at  .  .  . 

WSMR 

WSMR 

Balmertown,  Ontario,  Canada 

Transmitter  frequency, 

megahertz* 

2.2375 

2.6667 

2.6667 

Gyromagnetic  frequency, 

megahertz* 

1.404 

1.404 

1.638 

Propagation/f ield  line 

angle  <|),  degrees 

30.0 

30.0 

12.2 

*  Convert  to  radians  per 

second  for 

use  in  the 

computations . 
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The  polynomial  is  readily  differentiated  to  obtain  values  of  dy/dx  (at  line 
14),  and  division  by  G(x)  at  line  15  completes  the  determination  of  the  elec¬ 
tron  densities  over  the  considered  interval,  in  units  of  electrons  per  cubic 
centimeter. 

Results  for  the  example  run  are  plotted  in  Figure  6.4,  which  shows  altitude  in 
kilometers  on  the  vertical  axis,  and  the  logarithm  of  electron  density  (electrons 
per  cubic  centimeter)  on  the  horizontal  axis,  which  is  scaled  from  10*  to  10^. 
Numerical  values  of  electron  density  appear  to  the  right  of  the  altitude 
values . 


47 


SECTION  7 


An  Alternate  Procedure  for  Data  Processing 

The  data  reduction  procedure  based  on  Equation  (1.5),  and  implemented  in  the 
program  RLTIM3  described  in  Section  6,  is  what  we  shall  now  call  the  standard 
procedure  for  reducing  partial-reflection  data.  We  have  developed  another 
procedure  that  provides  benefits  not  otherwise  obtained,  and  this  non-standard 
procedure  is  the  subject  for  discussion  here. 

Several  investigators  (5,  6)  have  noted  a  strong  correlation  between  electron- 
density  profiles  N(h)  and  the  profile  of  ordinary  amplitude  AQ(h).  The  corre¬ 
lation  is  most  pronounced  at  low  altitude.  Belrose  and  Burke  (5)  used 
profiles  to  extrapolate  their  electron  density  profiles  to  altitudes  below 
which  results  were  otherwise  suspect.  They  drew  theoretical  support  for  this 
technique  from  a  paper  by  Booker  (7). 

Our  non-standard  method  is  designed  to  obtain  the  electron-density  profile 
more  or  less  directly  from  the  AQ(h)  profile.  First,  we  use  the  standard 
method  to  obtain  an  electron  density  profile  N(h)  that  is  restricted  to  mid¬ 
altitudes  (70-80  km  range)  only,  where  the  results  tend  to  have  maximum  valid¬ 
ity.  We  then  adjust  constants  and  in  this  expression, 

C, 

N'(h)  =  C2  •  h  •  AQ(h)  •  (£(h) )  \  (7.1) 

until  N'  agrees  as  closely  as  practical  with  N  over  the  altitude  range  of  N. 
Then  N*  is  taken  as  the  profile  of  electron  density  for  the  run. 

Kissick  (8)  has  developed  another  procedure,  more  elaborate  than  ours,  for 
obtaining  electron  densities  from  only  the  AQ(h)  profile.  His  work  lends 
strong  support  to  the  idea  that  the  ordinary  echoes  alone  are  rich  in  in¬ 
formation. 

The  rationale  for  Equation  (7.1)  will  now  be  explained. 
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The  function  E(h)  is  empirical,  and  is  given  in  Table  7.1.  It  has  unit  value 
at  low  altitude  and  increases  in  exponential  fashion  at  high  altitude.  The 
purpose  of  E(h)  is  to  compensate  for  the  fact  that  AQ(h)  is  reduced  at  high 
altitude  by  the  reflection  and  absorption  of  radio  waves  at  lower  altitude. 

Multiplication  of  by  altitude  h  in  Equation  (7.1)  is  done  for  this  reason: 

In  the  well-known  radar  equation,  echo  intensity  is  proportional  to  the  back- 

scattering  cross  section  of  the  target,  and  inversely  proportional  to  the 

4 

fourth  power  of  target  distance  (I  «  a/h  ).  The  cross  section  for  partial 

reflections  should  be  proportional  to  the  area  of  the  radar  beam,  which  varies 

2  2  2 
ash  (a  a  h  ) .  The  combining  of  these  relationships  gives  I  «  1/h  for  echo 

intensity,  and,  correspondingly,  Aq  «  1/h  for  echo  amplitude.  Multiplication 

of  Aq  by  h  in  Equation  (7.1)  is  thus  compensation  for  the  altitude  effect. 

Constants  and  in  Equation  (7.1)  are,  as  noted,  adjustable,  and  should  be 
selected  to  make  N’  follow  the  trend  of  N.  In  effect,  controls  the  ’’slope" 
of  N 1 ,  since  increasing  increases  the  value  of  N’  at  high  altitude  but  does 
not  change  the  value  at  low  altitude.  Constant  C^  simply  scales  the  overall 
result . 

An  example  from  the  non-standard  reduction  is  given  in  Figure  7.1.  The  Aq 
values  here  were  taken  from  column  3  of  Figure  6.2,  and  the  constants  were 
chosen  for  good  agreement  with  the  corresponding  standard  result  of  Figure  6.4. 
Constant  values  were  C^  =  3.0,  C^  =  0.14. 

What  are  the  advantages  of  the  non-standard  method?  In  the  first  place,  N’ (h) 
provides  electron  density  at  altitudes  lower  than  N(h)  can  reasonably  be 
extended,  which  is  the  same  extrapolation  of  results  used  by  Belrose  and 
Burke.  In  the  second  place,  N’ (h)  can  show  a  "fine  structure"  not  found  in 
N(h).  At  best,  the  standard  method  can  provide  only  gross  changes  in  electron 
density  with  altitude.  Any  interesting  features  that  cover  only  a  few  kilo¬ 
meters  will  be  lost.  However,  any  features  in  the  AQ(h)  profile  will  still  be 
present  in  the  N'(h)  profile. 

Results  are  given  in  Section  8.  In  all  cases,  results  were  obtained  using  the 
alternate,  non-standard  method. 
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TABLE  7.1.  TABULATIONS  OF  E(h) .  E(h)  IS  THE  EMPIRICAL  EXPONENTIAL  FUNCTION 
USED  IN  THE  NON-STANDARD  DATA  REDUCTION.  v(h)  IS  ELECTRON 
COLLISION  FREQUENCY  OBTAINED  FRCM  MEASUREMENTS  AT  WSMR  MADE  VIA 
THE  ROBIN  SPHERE  PROGRAM,  AND  USED  IN  THE  STANDARD  DATA  REDUCTION 


Altitude  h, 

km  ECh)  v(h) 


51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 


1.0000 
1.0001 
1.0001 
1.0002 
1.0003 
1.0005 
1.0007 
1.0010 
1.0015 
1.0020 
1.0028 
1.0037 
1.0049 
1.0065 
1.0085 
1.0109 
1.0139 
1.0175 
1.0218 
1.0268 
1.0328 
1.0396 
1.0474 
1.0564 
1.0664 
1.0777 
1.0904 
1.1045 
1 . 1202 
1.1377 
1.1570 
1.1785 
1.2023 
1.2288 
1.2582 
1.2910 
1.3276 
1.3684 
1 .4142 
1.4655 
1.5234 
1 .5887 
1.6627 
1.7468 
1 .8429 
1 .9531 
2.0800 
2.2270 
2.3982 
2.5987 


5.2000E7 

4.5000E7 

4.0000E7 

3.5000E7 

3. 1000E7 

2.6880E7 

2.3550E7 

2.0660E7 

1.8110E7 

1.5880E7 

1.4040E7 

1.2260E7 

1.0680E7 

9.3100E6 

8.  1050E6 

7.0460E6 

6.0890E6 

5.2270E6 

4.4940E6 

3.8700E6 

3.3410E6 

2 . 8860E6 

2.4870E6 

2.1380E6 

1.8320E6 

1.5600E6 

1.3230E6 

1. 1150E6 

9.4000E5 

7.8770E5 

6.6040E5 

5.5220E5 

4.6250E5 

3.8760E5 

3.2430E5 

2.7200E5 

2 . 2840E5 

1.9I80E5 

1.6160E5 

1.3640E5 

1.1500E5 

9.6930E4 

8. 1440E4 

6.7940E4 

5 . 6500E4 

4.6940E4 

3.8970E4 

3.3000E4 

2. 7000E4 

2.2500E4 
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Figure  7.  1.  Example  result,  non-standard  data  processing. 


SECTION  8 


i 


Results 

From  this  partial  reflection  experiment,  682  electron  density  profiles  have 
been  produced.  Because  this  number  is  so  large,  the  profiles  are  contained  in 
a  separate  volume.  Figure  8.1  shows  the  format  we  have  used  to  present  the 
profiles  determined  from  the  individual  runs.  In  each  case  there  is  a  tabula¬ 
tion  of  electron  density  at  1-kilometer  increments,  and  a  semilogarithmic 
plot  of  electron  density  vs.  altitude.  Recall  that  the  time  duration  of  a  run 
is  on  the  order  of  10  minutes,  ar.d  that  results  are  based  on  a  time  average 
of  the  several  hundred  measurements  of  echo  amplitude  made  during  the  run. 

Tables  8.1,  8.2,  and  8.3  list  the  date  and  start-times  of  the  runs. 

The  41  runs  in  Table  8.1  were  made  at  WSMR  in  1975  by  Mr.  Glenn  D.  Falcon  of 
the  Institute  for  Telecommunication  Sciences,  Office  of  Telecommunications, 

U.  S.  Department  of  Commerce,  Boulder,  Colorado.  These  data  were  taken  prior 
to  PSL  participation  in  the  experiment.  However,  data  reduction  for  these 
runs  was  performed  by  PSL. 

Table  8.2  lists  413  runs  made  by  PSL  at  WSMR  between  September  1977  and  August 
1979. 

Table  8.3  lists  198  runs  made  by  PSL  at  Balmertown,  Ontario,  Canada,  in  February 
1979.  The  station  was  moved  to  Canada  as  part  of  the  Atmospheric  Sciences 
Laboratory  scientific  program  for  the  solar  eclipse  of  26  February  1979. 

Data  from  the  eclipse  program  spanning  totality  were  also  processed  in  30 
1-minute  segments,  extending  from  1045  to  1115  local  time*  (Totality  was 
from  approximately  1053  to  1055.)  These  30  profiles  are  included  in  the 
results  volume. 

Figures  8.2  -  8.11  contain  averages  of  results  for  the  WSMR  data  on  a  monthly 
basis.  In  forming  the  averages,  we  used  only  data  recorded  during  the  midday 
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Figure  8.1.  Key  to  electron  density  profiles. 


hours  of  1000  to  1400  MST,  and  only  covered  those  months  for  which  a  signifi¬ 
cant  number  of  midday  runs  were  available.  Within  the  plots,  'o’  marks  the 
average  value,  while  the  marks  are  separated  from  the  average  by  one 
standard  deviation. 

In  regard  to  the  eclipse  data,  solar  photons  are  the  primary  agent  for  ioniza¬ 
tion  of  the  atmosphere  in  the  D  region,  and  a  total  solar  eclipse  can  allow 
measurements  showing  the  dynamics  of  the  ionization  process.  However,  the  sun 
was  very  active  during  February  1979,  and,  apparently,  particles  in  the  solar 
wind  were  a  major  factor  in  D-region  ionization.  Thus  changes  in  ionization 
during  the  eclipse  were  not  pronounced.  Figure  8.12  shows  time  variation  in 
electron  density  obtained  from  partial  reflections  during  the  eclipse,  at 
altitudes  of  60,  65,  70,  75  km.  A  small  decrease  of  electron  density  at  the 
lower  two  altitudes  is  seen,  but,  at  the  higher  two  altitudes,  there  appears 
to  be  an  increase  in  electron  density.  Interpretation  of  these  results  is 
still  pending. 
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TABLE 


Time,  Solar 
Mtn.  Zenith 

Date  Daylight  Angle 

0  9- 26 -"7  5  05  57  - 

0622  98. 6  * 

0633  96.3 

0649  92.9 

0659  90.8 

0718  86.8 

0735  83.2 

0747  80.7 

0810  75.9 

0836  70.6 

10-01-75  1330  36.3 

1355  37.9 

1415  39.8 

1445  43.4 

1520  48.7 

1545  52.9 

1635  62.1 

1713  69.6 

1800  79.2 

1840  87.5 

1955 

10-02-75  1205  38.3 

1231  36.6 

1311  36.1 

1410  39.6 

1435  42.5 

1455  45.2 

1516  48.3 

1531  50.8 

1550  54.0 

1613  58.2 

1700  67.2 

1800  79.4 

1900  92.0 

2010 

10-03-75  1040  49.2 

1140  41.0 

1210  38.2 

1237  36.8 

1412  40.2 

1436  42.9 


.1.  RUN  INDEX  FOR  WSMR,  1975 


Sequence 

No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


TABLE  8.2, 

RUN  INDEX  FOR 

WSMR,  1977- 

-1979 

Solar 

Local 

Zenith 

Sequence 

Date 

Time*  Angle 

No. 

09-20-77 

1135 

37.  3  ° 

1 

10-12-77 

1230 

40,6° 

41 

09-23-77 

1135 

38.3 

2 

1300 

40.0 

42 

1205 

35.  1 

3 

1315 

40.  2 

43 

1713 

67,  8 

4 

1331 

40.7 

44 

1735 

72.  4 

5 

1400 

42.  5 

45 

09-26-77 

0830 

71.  9 

6 

1426 

45.  0 

46 

0930 

59.  9 

7 

1443 

46.  9 

47 

1030 

49.  0 

8 

1505 

49.  9 

48 

1130 

39.  9 

9 

1525 

52.  9 

49 

1230 

34.6 

10 

1549 

56.  8 

50 

1330 

34.6 

11 

10-  19-77 

0930 

65.  7 

51 

1430 

40.  0 

12 

1030 

55.6 

52 

1515 

46.5 

13 

1130 

47.  7 

53 

09-27-77 

1015 

51.8 

14 

1230 

43.2 

54 

1253 

3  4.3 

15 

1330 

43.  2 

55 

09-28-77 

1253 

34.6 

16 

1430 

47.  7 

56 

09-29-77 

0755 

79.8 

17 

10-25-77 

1253 

44.  7 

57 

0845 

60.6 

18 

10-26-77 

0953 

63.4 

58 

0953 

56.4 

19 

1053 

54.  3 

59 

1053 

46.  i 

20 

1153 

47.  8 

60 

1153 

38.4 

21 

1253 

45.  0 

61 

1253 

3  5.0 

22 

1353 

46.  8 

62 

1400 

37.8 

23 

170  5 

73.  8 

63 

1415 

3  9.3 

24 

10-31-77 

1155 

46.  7 

64 

1430 

4  1.0 

25 

11-01-77 

1153 

47.  0 

65 

1445 

43.0 

26 

11-02-77 

0855 

64.  7 

66 

1500 

45.  1 

27 

0930 

59.  3 

67 

1535 

50.7 

28 

1000 

55.  3 

68 

1551 

53.  5 

29 

1030 

51.9 

69 

1609 

56.8 

30 

1100 

49.4 

70 

09-30-77 

1255 

35.4 

31 

1130 

47.  8 

71 

10-05-77 

1253 

37.  4 

32 

1200 

47.3 

72 

o 

t 

o 

( 

1253 

39.  3 

33 

1230 

47.  8 

73 

10-11-77 

1253 

39.  7 

34 

1300 

49.4 

74 

1320 

39.  9 

35 

1330 

52.0 

7  5 

o 

i 

i 

0833 

74.  9 

36 

1403 

55.  7 

76 

0930 

63.  9 

37 

1430 

59.4 

77 

1115 

47.2 

38 

1500 

64.0 

78 

1130 

45.4 

39 

11-03-77 

0900 

64.  1 

79 

1215 

41,4 

40 

0930 

59.6 

80 

*  Mountain  daylight  time  thru  10-26-77, 

Mountain  standard  time  from  10-31-77. 

Note:  Mountain  daylight  times  on  the  data  plots  are  marked  with  a  . 
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r - 

Table  8.2  (cont) 

Solar 

Local 

Zenith 

Sequence 

Date  Time* 

Angle 

No. 

11-03-77  1001 

55.4° 

81 

11-05-77 

0615 

95.  3° 

121 

1030 

52.2 

82 

0645 

89.2 

122 

1100 

49.  7 

83 

0715 

83.  3 

123 

1131 

48.  1 

84 

0745 

77.6 

124 

1200 

47.6 

85 

0815 

72.2 

125 

1230 

48.  1 

86 

0831 

69.4 

126 

1300 

49.  7 

87 

0900 

64.  6 

127 

1330 

52.  2 

88 

0930 

60.  1 

128 

1400 

55.6 

89 

1000 

56.  1 

129 

1430 

59.6 

90 

1030 

52.  8 

130 

1500 

64.  2 

91 

1100 

50.  3 

131 

1530 

69.  2 

92 

1130 

48.  7 

132 

11-04-77  1225 

48.  3 

93 

1200 

48.2 

133 

1300 

50.0 

94 

1230 

48.  7 

134 

1330 

52.  5 

95 

1300 

50.  3 

135 

1400 

55.  9 

96 

1330 

52.  8 

136 

1430 

59.  9 

97 

1400 

56.  1 

137 

1500 

64.4 

98 

11-09-77 

0930 

61.0 

138 

1530 

69.4 

99 

1000 

57.  1 

139 

1600 

74.  7 

100 

1030 

53.8 

140 

1630 

80.  3 

101 

1100 

51.4 

141 

1645 

83.  2 

102 

1130 

49.9 

142 

1700 

86.  2 

103 

1200 

49.  3 

143 

1715 

89.  1 

104 

1230 

49.9 

144 

1730 

92.  1 

105 

1300 

51.4 

145 

1745 

95.  2 

106 

1330 

53.  9 

146 

1830 

- 

107 

1400 

57.  1 

147 

:  1900 

- 

108 

1430 

61.0 

148 

2000 

- 

109 

1500 

65.5 

149 

2100 

- 

110 

11-22-77 

0945 

61.7 

150 

2200 

- 

111 

11-23-77 

1130 

53.2 

151 

230u 

- 

112 

1200 

52.  8 

152 

11-05-77  0000 

- 

113 

1230 

53.  3 

153 

0100 

- 

114 

1300 

54.  7 

154 

0200 

115 

1345 

58.  5 

155 

0300 

- 

116 

1415 

61.9 

156 

0400 

- 

117 

1445 

65.  9 

157 

0500 

- 

118 

1520 

71.2 

158 

0545 

- 

119 

12-07-77 

0930 

65.6 

159 

0600 

98.  3 

120 

1000 

62.0 

160 

*  Mountain 

standard 

time. 
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Table  8.2  (cont) 


Solar 

Local  Zenith  Sequence 

Date  Time*  Angle  No. 

12-07-77  1030  59.  1*  16 1 

1100  56.8  162 

1130  55.5  163 

1200  55.0  164 

1230  55.5  165 

1300  56.9  166 

1330  59. 1  167 

1400  62.0  168 

1430  65.  7  169 

1500  69.8  170 

12-14-77  1015  61.0  171 

1045  58.4  172 

1115  56.7  173 

1145  55.7  174 

1215  55.  7  175 

1245  56.7  176 

1300  57.4  177 

1315  58.4  178 

1330  59.6  179 

1400  62.6  180 

12-21-77  1015  61.3  181 

1044  58.8  182 

1115  56.9  183 

1145  56.0  184 

1215  56.0  185 

1245  56.9  186 

1315  58.7  187 

1345  61.  3  188 

1415  64.5  189 

1445  68.4  190 

12-28-77  1130  56.3  191 

1200  55.8  192 

1230  56.3  193 

01-04-78  0915  67.9  194 

0945  64. 1  195 

1015  60.8  196 

1045  58.2  197 

1115  56.4  198 

1145  55.5  199 

1215  55.5  200 


01-04-78  1245  56.4  201 

1315  58.2  202 

1345  60.8  203 

1415  64.0  204 

1445  67.  9  205 

01-11-78  1015  60.1  206 

1045  57.4  207 

1115  55.6  208 

1145  54.7  209 

1215  54.  7  210 

1345  60.0  211 

1415  63.4  212 

1445  67.3  213 

01-18-78  0935  63.7  214 

1000  60.6  215 

1100  55.3  216 

1130  53.9  217 

1230  53.  9  218 

1300  55.3  219 

1330  57.6  220 

1400  60.6  221 

1430  64.3  222 

1500  68.6  223 

01-25-78  0930  63.2  224 

1000  59.4  225 

1030  56.3  226 

1100  53.9  227 

1130  52.5  228 

1200  52.0  229 

1230  52.5  230 

1300  53.9  231 

1330  56.2  232 

1400  59.4  233 

1430  63. 1  234 

02-01-78  0900  66.2  235 

0930  61.8  236 

1000  57.9  237 

1030  54.7  238 

1100  52.3  239 

1130  50.8  240 


*  Mountain  standard  time. 
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Table  8.2  (cont) 


Date 

02-01 


02-08 


02-15 


02-22 


03-01 


03-15 


03-22 


Solar 

Local  Zenith 

Time*  Angle 

78  1200  50.  3° 

1230  50.8 

1300  52.3 

1330  54.7 

1400  57.9 

1430  61.7 

1520  69.3 

78  0920  61.6 

1000  56.2 

1030  52.9 

1100  50.4 

1345  54.4 

1415  58.1 

1445  62.3 

78  0900  63.  1 

1000  54. 3 

1030  50.9 

1100  48.3 

1225  46.5 

1335  51.4 

1415  56.3 

1445  60.6 

78  1045  47.3 

1115  45.1 

1145  43.9 

1215  43.9 

1245  45.0 

78  1240  42.3 

1310  44.5 

1330  46.5 

1355  49.5 

78  1210  36.0 

1230  36.6 

1300  38.6 

1330  41.8 

1400  45.8 

78  1045  37.6 

1115  34.9 

1145  33.4 

1215  33.4 


Sequence 

No. 

241  03-22-78 

242 

243 

244 

245 

246  03-29-78 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256  04-05-78 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266  04-13-78 

267 

268 

269 

270 

271  04-14-78 

272 

273 

274 

275 

276 

2  n 

278 

279 

280 


1245 

. 

00 

o 

281 

1315 

37.6 

282 

1345 

41.4 

283 

1415 

46.0 

284 

1445 

51.  1 

285 

1000 

41.5 

286 

1030 

37.  1 

287 

1100 

33.6 

288 

1130 

31.2 

289 

1200 

30.4 

290 

1230 

31.2 

291 

1300 

33.  5 

292 

1330 

37.0 

293 

1400 

41.5 

2  94 

1430 

46.  5 

295 

1000 

39.4 

2  96 

1030 

34.  8 

297 

1100 

31.  1 

298 

1130 

28.6 

299 

1200 

27.  7 

300 

1230 

28.6 

301 

1300 

31.0 

302 

1330 

34.  8 

303 

1400 

39.4 

304 

1430 

44.6 

305 

1445 

45.  5 

306 

1515 

51.4 

307 

1545 

57.  5 

308 

1615 

63.  7 

309 

1645 

70.0 

310 

1043 

30.  2 

311 

1115 

26.  5 

312 

1145 

24.  5 

313 

1215 

24.  5 

314 

1245 

26.4 

315 

1315 

29.  9 

316 

1345 

34.4 

317 

1415 

39.6 

318 

1445 

45.2 

319 

1515 

51.2 

320 

*  Mountain  standard  time. 
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Table  8.2  (cont) 


Solar 


Local 

Zenith 

Sequence 

Date 

Time  * 

Angie 

No. 

07- 19-78 

1315 

11.8° 

321 

08-18-78 

0955 

46.  3° 

361 

1400 

17.  5 

322 

1025 

40.2 

362 

1430 

23.0 

323 

1055 

34.  3 

363 

1500 

28.  9 

324 

1125 

28.  9 

364 

07-21-78 

1445 

26.  1 

325 

1155 

24.2 

365 

08-07-78 

1100 

31.  3 

326 

1225 

20.  7 

366 

1130 

25.  7 

327 

1255 

19.  2 

367 

1200 

21.2 

328 

1325 

20.  0 

368 

1250 

15.  9 

329 

1400 

23.6 

369 

1310 

15.9 

330 

12-21-78 

1  125 

56.  5 

370 

o 

00 

1 

4^ 

i 

00 

1030 

38.  5 

331 

1150 

55.  9 

371 

1100 

32.6 

332 

01-22-79 

1200 

52.  7 

372 

1200 

22.  5 

333 

1215 

52.8 

373 

1230 

20.0 

334 

1230 

53.2 

374 

1300 

17.  9 

335 

05-25-79 

1130 

23.4 

375 

1330 

19.  1 

336 

1200 

18.0 

376 

1400 

22.  5 

337 

1230 

13.8 

377 

1430 

27.  2 

338 

1300 

12.0 

378 

08- 16-78 

0900 

57.4 

339 

1330 

13.  7 

379 

0930 

51.  1 

340 

06-29-79 

1000 

40.  5 

380 

1000 

44.  9 

341 

1030 

34.  2 

381 

1030 

38.  8 

342 

1100 

27.  9 

382 

1100 

33.0 

343 

1130 

21.8 

383 

1130 

27.6 

344 

1200 

16.  0 

384 

1200 

23.  0 

345 

1230 

11.2 

385 

1230 

19.  7 

346 

1300 

9.0 

386 

1300 

18.  5 

347 

1330 

11.2 

387 

1330 

19.7 

348 

1400 

16.0 

388 

1400 

23.0 

349 

07-09-79 

1150 

18.  3 

389 

08- 17-78 

0900 

57.  5 

350 

1230 

11.8 

390 

0925 

52.  3 

351 

1320 

11.8 

391 

0950 

47.  1 

352 

1400 

16.  5 

392 

1025 

40.  0 

353 

1424 

2  1.0 

393 

1055 

34.  1 

3  54 

1445 

25.  2 

3  94 

1125 

28.  7 

355 

07- 10-79 

1120 

24.2 

395 

1155 

24.0 

356 

1153 

17.  8 

3  96 

1225 

20.  4 

357 

1240 

10.  8 

397 

1250 

19.0 

358 

1400 

16.6 

398 

1330 

20.  0 

359 

1430 

22.  2 

399 

08-18-78 

0935 

50. 4 

360 

07-11-79 

1000 

40.9 

400 

*  All  times  are  mountain  daylight,  except 
on  12-21-78  and  01-22-79,  which  are 
mountain  standard. 
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Table  8*2  (cont) 


Date 

07-11-79 


08-01-79 


*  Mountain 


Local 

Solar 

Zenith 

Sequence 

Time* 

Angle 

No. 

1030 

34.  6  0 

401 

1100 

28.  3 

402 

1130 

22.  3 

403 

1200 

16.6 

404 

1230 

12.0 

405 

1315 

10.  5 

406 

1345 

14.  1 

40  7 

1430 

22.  3 

408 

1625 

47.  9 

409 

1700 

55.  3 

410 

1740 

63.  7 

411 

1815 

71.  1 

412 

1845 

77.  3 

413 

daylight  time 
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TABLE  8.3.  RUN  INDEX  FOR  BALMERTOWN,  CANADA,  1979 


Time,  Solar 

Central  Zenith  Sequence 


Date 

Std. 

Angle 

No 

0*3-08-79 

1100 

68.  3 

1 

1155 

67.  1 

2 

1 320 

69.  3 

3 

1355 

71.6 

4 

1455 

77.  1 

5 

1645 

91.  1 

6 

02-09-79 

0850 

78.6 

7 

1130 

67.  1 

8 

1700 

93.0 

9 

02-10-79 

1305 

67,  9 

10 

1355 

71.0 

11 

1510 

78.  3 

12 

02-11-79 

0930 

73.  8 

13 

1000 

71.2 

14 

1050 

67.  9 

15 

1200 

66.  1 

16 

1300 

67.  4 

17 

02-  12-79 

0800 

84.0 

18 

0900 

76.  7 

19 

22  15 

- 

20 

02-13-79 

0820 

81.  1 

21 

0900 

76.  4 

22 

1000 

70.6 

23 

1050 

67.  3 

24 

1150 

65.6 

25 

1335 

68.  7 

26 

1400 

70.  5 

27 

1500 

76.  3 

28 

1520 

78.6 

29 

1600 

83.6 

30 

02-16-79 

1050 

66.  3 

31 

1150 

64.6 

32 

1210 

64.6 

33 

1250 

65.  5 

34 

1350 

68.  8 

35 

1450 

74.  3 

36 

1600 

82.  8 

37 

1700 

91.3 

38 

1720 

94.  3 

39 

1810 

- 

40 

02-17-79 

0850 

76.4 

41 

0910 

74.  1 

42 

0950 

70.2 

43 

1050 

66.  0 

44 

1150 

64.  3 

45 

1350 

68.  5 

46 

1410 

70.2 

47 

1450 

74.  1 

48 

1550 

81,2 

49 

1650 

89.  5 

50 

02-17-79  1730 

95.  5 

51 

1740 

97.0 

52 

1815 

- 

53 

02-18-79  0845 

76.  7 

54 

0905 

74.4 

55 

0950 

69.  9 

56 

1050 

65.7 

57 

1150 

63.  9 

58 

1250 

64.  8 

59 

1350 

68.2 

60 

1425 

71.2 

61 

02-19-79  0830 

78.  2 

62 

0905 

74.  1 

63 

0950 

69.6 

64 

1050 

65.4 

65 

1150 

63.6 

66 

1250 

64.  5 

67 

1350 

67.  9 

68 

1415 

70.0 

69 

1450 

73.  5 

70 

1550 

80.6 

71 

1650 

89.0 

72 

02-20-79  0630 

94.  8 

73 

0700 

90.  3 

74 

0730 

85.  9 

75 

0850 

75.  5 

76 

0910 

73.2 

7? 

0950 

69.  3 

78 

1050 

65.  0 

79 

1150 

63.3 

80 

1250 

64.  1 

81 

1350 

67.6 

82 

1450 

73.  1 

83 

1550 

80.  4 

84 

1645 

88.  0 

85 

02-21-79  0630 

94.6 

86 

0700 

90.0 

87 

0730 

85.  7 

88 

0850 

75.2 

89 

0  950 

68.  9 

90 

1010 

67.  3 

91 

1050 

64.  7 

92 

1150 

62. 9 

93 

1250 

63.  8 

94 

1350 

67.2 

95 

1450 

72.  8 

96 

1550 

80.  1 

97 

32-22- 79  0530 

- 

98 

0600  98. 9  99 

0630  94.3  100 
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Table  8.3  (cont) 


Date 

02-22-79 


02-23-79 


02-24-79 


02-25-79 


Time, 

Solar 

Central 

Zenith 

Sequence 

Std. 

Angle 

No* 

0700 

89.  8 

101 

02-25-79 

1115 

62.2 

151 

0730 

85.  4 

102 

1130 

61.8 

152 

0850 

74.  9 

103 

1145 

61.6 

153 

0950 

68.  6 

104 

1200 

61.5 

154 

1020 

66.  2 

105 

1215 

61.5 

155 

1050 

64.  3 

106 

1300 

62.  8 

1 56 

1130 

62.6 

107 

02-26-79 

0540 

- 

157 

1250 

63.  4 

108 

0600 

97.  9 

158 

1350 

66.  9 

109 

0630 

93.2 

159 

0530 

1 10 

0645 

90.  9 

160 

0600 

98.  7 

111 

0700 

88.6 

161 

0630 

94.  0 

112 

0716 

86.  3 

162 

0700 

89.  5 

113 

0800 

80.0 

163 

0730 

85.  1 

1  14 

0830 

76.  1 

164 

0850 

74.  6 

115 
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Figure  8.2.  Mean  and  standard  deviation  electra  density  profile  for  September  1977. 
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Figure  8.5.  Mean  and  standard  deviation  electra  density  profile  for  December  1977. 
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Figure  8.8.  Mean  and  standard  deviation  electra  density  profile  for  March  1977. 


Figure  8.9.  Mean  and  standard  deviation  electra  density  profile  for  April  1977. 
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Mean  and  standard  deviation  electra  density  profile  for  May  1977. 
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APPENDIX  I 


Schematic  Diagrams 

The  following  schematics  for  the  partial  reflection  sounder  are  contained 
in  this  Appendix: 
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APPENDIX  II 


Averaging  Program 

This  appendix  contains  a  listing  of  the  FORTRAN  program  used  to  time-average 
the  data.  The  program  is  run  at  White  Sands  Missile  Range  on  the  UNIVAC  1108 
system. 

Pages  2-5  contain  the  main  program,  named  FLOOD.  Subroutines  INDATA3  and 
UNPACK  are  found  on  pages  6  and  7.  Page  8  contains  information  regarding 
input  cards . 
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START  Of  LOOP  TOR  DATA  ACCUMULATION 
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-  cjitQtu _ kiiQUftisift- 


SUfiROUTlNC  UNPACK  <  JOAT  A  ,  |NotO  ,*4UH»  ooooo  t  Oo 

0  1  MpHS  I  ON  J0ATAU1  .fMOLDIU  00000200 

,C*n  00000300 

o  o  i  vxr  jtttru^  mnwiirron 


I.  Tape  control.  One  card,  format  215,3011 


Column 

Format 

Name 

Description 

1-5 

15 

NOFILE 

Number  of  files  (N)  to  skip  on  the  front 
end  of  the  data  tape.  Note:  the  first 
file  contains  no  data,  so  ,ffile  1"  of  data 
is  actually  the  second  file  on  the  tape. 

6-10 

15 

KFILE 

Number  of  files  (M) ,  after  the  N  skipped 
files,  to  be  considered  for  processing 
(including  any  files  not  to  be  processed 
via  flags  in  the  following  M  columns). 

11-40 

3011 

IFILE 

One  digit  for  each  of  the  M  files  to  be 
considered;  column  11  for  the  first,  column 
10+M  for  the  last.  Files  flagged  with  a  0 
(or  blank)  will  be  processed,  those  with  a 

1  will  be  skipped  over. 

II.  Receiver  calibration  array  RX;  64  values  on  8  cards,  8  numbers  per  card, 
each  of  format  8(3X,F7.4). 

III.  Information  relating  to  the  files  to  be  processed.  One  card  per  file  of 
format  415 , 7A6 , 13, 215 , 12. 


Column 

Format 

Name 

Description 

1-5 

15 

JFILE 

File  number  to  be  printed  on  the  listing. 

6-10 

15 

ISTHT 

Starting  height  of  echo  samples. 

11-15 

15 

NSEG 

Number  of  records  per  segment. 

16-20 

15 

IJKHT 

Reference  height  index  for  noise  thresholds 

21-62 

7A6 

I DATE 1 

through  IDATE7 

Date,  time,  and  comments  to  be  printed  on 
listing.  Date  and  time  are  also  used  to 
compute  solar  zenith  angle. 

63-65 

13 

IPLT 

•  Plot  flag;  0  for  plots,  1  for  no  plots  of 
average  amplitudes. 

66-70 

15 

MAXI 

First  noise  threshold  for  screening  raw 
data . 

71-75 

15 

MAX2 

Second  noise  threshold  for  screening  raw 
data. 

76-77 

12 

I  PUN 

•  Punch  flag;  0  for  cards,  1  for  no  cards 
containing  the  average  amplitudes. 

•  Note:  we  generally  do  not  call  for  cards 
or  plots . 
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Solar  Zenith  Angle 

Solar  zenith  angle,  which  is  the  angular  distance  of  the  sun  away  from  the 
local  vertical,  is  of  importance  in  this  project  because  of  the  fact  that 
changes  in  this  angle  are  generally  well  correlated  with  changes  in  electron 
density  in  the  D  region. 

We  describe  here  the  simple  procedure  we  have  used  to  estimate  solar  zenith 
angle  as  a  function  of  time  and  date.  (Accuracy  of  the  estimate  is  sufficient 
for  our  purpose.)  Parameters  required  are  the  latitude  of  the  station  site, 
and  a  factor  based  on  the  time  of  the  most  recent  autumnal  equinox.  The 
computer  program  ZEN  in  Figure  2,  which  is  written  in  APL,  performs  the  required 
computations.  Comments  follow  that  are  keyed  by  number  to  the  line  numbers  in 
Figure  2. 


The  basic  calculation  involves  the  spherical  triangle  on  the  celestial  sphere 
shown  in  Figure  1. 


Figure  1.  Solar  Zenith  Angle  Triangle. 


Verticies  N,  V,  and  S  of  the  triangle  correspond  respectively  to  projections 
of  the  earth’s  spin  axis,  local  vertical  at  the  station  site,  and  the  line 
from  earth  to  sun.  The  dotted  line  is  the  projection  of  the  earth’s  equitorial 
plane. 
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Angle  LAT  (between  points  A  and  V)  is  the  station  latitude,  LHA  (between  sides 
NV  and  NS)  is  the  local  hour  angle,  DEC  (between  points  S  and  B)  is  the  solar 
declination,  and  SZA  (between  points  V  and  S)  is  the  desired  solar  zenith 
angle. 

Angle  LHA  grows  uniformly  with  time  at  the  rate  of  360°  per  day,  and  is  zero 
at  local  noon.  Small  uncertainties  exist  here  because  of  the  fact  that  LHA  = 

0  does  not  generally  occur  at  precisely  noon  local  time,  and  may  occur  as  much 
as  10  or  15  minutes  before  or  after  noon.  However,  the  corresponding  un¬ 
certainty  in  solar  zenith  angle  is,  for  our  purpose,  negligible.  Note  that 
one  should  always  employ  local  standard  time  and  not  daylight-savings  time! 

Angle  DEC  is  a  sinusoidal  function  of  time,  with  an  amplitude  of  23.5°  and  a 
period  of  one  year.  It  is  zero  at  the  equinoxes,  and  changes  from  negative  to 
positive  at  the  autumnal  equinox.  (Figure  1  shows  a  positive  DEC.) 

Angle  SZA  is  found  using  this  standard  equation  of  spherical  trigonometry: 

cos  SZA  =  sin  LAT  •  sin  DEC  +  cos  LAT  *  cos  DEC  •  cos  LHA.  (1) 

Program  comments  (for  Figure  2) 

6.  As  described  in  the  comments  of  line  1-5,  BASE  is  to  be  set  to  the 
time  interval  in  hours  between  the  autumnal  equinox  preceding  the  time  of 
interest  and  00:00  of  the  subsequent  1  October  (all  times  local  standard 
time).  The  value  of  BASE  given  here  is  for  the  autumnal  equinox  of  1977  which 
occurred  at  20:30  mountain  standard  time  on  22  September  1977.  From  this 
event  to  00:00  hours  of  23  September  is  3.50  hours,  and  from  00:00  of  23 
September  to  00:00  hours  of  1  October  is  24  •  8  =  192  hours:  Thus,  BASE  = 

3.50  +  192  =  195.5  hours.  (See  Table  1  for  other  BASE  values.) 

7.  Our  station,  at  WSMR,  is  located  at  32°  24’  north  latitude.  (At  the 
Ontario,  Canada  site,  the  latitude  was  51°  5'.) 

8.  Days  per  month  from  October  1977  through  September  1978,  inclusive. 
The  time  of  interest  must  occur  in  one  of  these  months. 
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9.  Conversion  of  latitude  in  degrees,  minutes,  and  seconds  to  radians. 

10.-13.  Month,  day,  hour,  and  minute  of  the  time  of  interest.  Note  that 
these  4  values  are  entered  into  the  argument  of  the  call-statement  for  program 
ZEN. 

14.  This  statement  converts  the  month  value  from  the  calendar  scale  in 
which  1  =  January  to  the  "fiscal”  scale  in  which  1  =  October. 

15.  Month  is  converted  into  total  elapsed  days;  0  for  October,  31  for 
November,  61  for  December,  92  for  January,  etc. 

16.  Using  the  result  from  line  15,  the  time  of  interest  is  converted 
into  the  number  of  hours  (At)  from  the  autumnal  equinox  to  the  time  of  in¬ 
terest  . 

17.  SUN  is  an  angle  proportional  to  At: 

SUN  =  2n  •  (At  hours)  t  (number  of  hours  in  one  year)  radians 

18.  Local  hour  angle  LHA: 

LHA  =  2n  •  (hour  of  day  -  12  hours)  t  24  hours  radians 

19.  The  solar  declination  (DEC)  is  23.5°  times  the  sine  of  angle  SUN 
(from  line  17).  SDEC  is  the  sine  of  DEC. 

20.  CDEC  is  the  cosine  of  DEC. 

21.  Solar  zenith  angle  computed  via  Eq.  (1)  above. 

22.  Printout  of  solar  zenith  angle. 
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Table  1 .  Base  values  for  zenith  angle  calculations 


Values  applicable  for  WSMR;  mountain  standard  time 


from 

to 

base 

Oct  1974 

Sept  1975 

189.00  hours 

Oct  1975 

Sept  1976 

183.05 

Oct  1976 

Sept  1977 

201.18 

Oct  1977 

Sept  1978 

195.50 

Oct  1978 

Sept  1979 

189.57 

Oct  1979 

Sept  1980 

183.71 
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APPENDIX  IV 


Receiver  Calibration 


The  receiver  used  in  the  partial  reflection  sounder  is  calibrated  by  applying 
the  cw  signal  from  an  rf  signal  generator  to  the  receiver  input  and  observing 
the  digitized  values  of  output  signal  strength  vs.  input  amplitude  (see  Section 
5.5).  The  receiver  should  be  calibrated  following  any  adjustment  that  could 
influ. nee  its  sensitivity,  or  after  any  period  of  time  (weeks,  typically)  in 
which  the  sensitivity  might  drift. 

In  practice,  digitized  values  obtained  in  the  calibration  show  a  small  random 
fluctuation  with  time,  and  so  the  average  of  many  values  is  used  in  each  case. 
With  constant  input  to  the  receiver,  the  digitized  signal  strength  is  recorded 
on  magnetic  tape  for  several  minutes  in  the  same  format  as  reflection  data, 
and  the  program  used  to  average  the  reflection  data  is  also  used  to  average 
the  calibration  data.  A  "run"  of  calibration  data  is  identified  with  a  par¬ 
ticular  input  level. 

Table  1  shows  typical  results  of  input  signal  amplitude  (in  microvolts)  vs. 
average  counts.  Each  change  in  input  voltage  corresponded  to  a  1  dB  change  in 
input  power.  The  dBm  input  values  are  included  in  the  table  only  for  reference, 
and  play  no  part  in  the  calibration  procedure. 

We  assume  that  input  amplitudes  from  the  signal  generator  can  be  expressed 
as  a  function  of  average  counts  C^,  with  the  relationship  being  of  the  polyno¬ 
mial  form 


S. 

i 


3 

1  A.(C.)J, 

j=0  J  1 


(1) 


where  A^,  A^,  and  A^  are  constants.  Best  values  of  the  (over  determined) 

constants  are  found  using  the  method  of  least  squares.  For  the  values  in 
Table  1,  best  values  are  as  follows: 
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TABLE  1.  TYPICAL  RECIEVER  CALIBRATION  DATA 


Power  input 

to  receiver  Average 


from  signal 

generator, 

dBm 

Input 

amplitude  (S), 
microvolts* 

digitized  count 
at  receiver 
output 

-97 

3.158 

57.97 

-98 

2.815 

48.06 

-99 

2.509 

39.51 

-100 

2.236 

32.43 

-101 

1.993 

26.23 

-102 

1.776 

24.05 

-103 

1.583 

19.63 

-104 

1.411 

16.13 

-105 

1.257 

13.74 

-106 

1.121 

11.35 

-107 

0.998 

9.63 

-108 

0.890 

8.09 

-109 

0.793 

7.05 

-110 

0.707 

6.03 

-111 

0.630 

5.14 

-112 

0.562 

4.10 

-113 

0.500 

3.38 

-114 

0.446 

2.72 

-115 

0.398 

2.09 

-116 

0.354 

1.54 

-117 

0.316 

0.99 

Input  amplitude  is  computed  from  input  power,  and  not  measured. 
It  is  the  input  power  that  is  measured. 
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AQ  =  0.21990 
A2  =  0.088315 
A2  =  -0.0010286 
A3  =  6.5537  x  10' 


When  it  comes  to  reducing  the  partial-reflection  data,  we  require  the  ratio  of 
received  amplitudes;  extraordinary-mode  amplitude  divided  by  ordinary-mode 
amplitude.  It  is  by  means  of  the  receiver  calibration  "curve"  that  average 
digital  counts,  which  constitute  the  data,  are  converted  into  corresponding 
signal  amplitude.  Any  linear  factor  introduced  into  the  calibration  will 
cancel  in  the  ratio.  Thus,  and  only  for  convenience,  we  scale  the  constants 
so  that  a  maximum  count  of  63  corresponds  to  an  amplitude  of  63  (units  arbi¬ 
trary).  After  each  constant  is  multiplied  by  the  appropriate  scale  factor, 
these  values  obtain: 

Aq  =  4.1480 
Aj  =  1.6659 
A2  =  -0.019404 
A3  =  1.2362  x  It)"4. 

If  we  use  these  constants  in  Equation  (1),  signal  amplitude  (S)  as  a  function 
of  count  (C)  is  as  shown  in  Table  2.  To  use  this  calibration,  the  64  values 
in  the  "S"-columns  of  Table  2  are  entered  into  the  data  reduction  program 
(Appendix  II)  as  the  array  RX. 
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TABLE  2.  TABULATED  RECEIVER  CALIBRATION 
DATA  OF  TABLE  1 


'CURVE"  OBTAINED 


ount  (C) 

Amplitude  (S) 

(c: 

0 

0.0000 

32 

1 

5.7947 

33 

2 

7.4032 

34 

3 

8.9745 

35 

4 

10.5091 

36 

5 

12.0080 

37 

6 

13.4717 

38 

7 

14.9011 

39 

8 

16.2968 

40 

9 

17.6597 

41 

10 

18.9905 

42 

11 

20.2898 

43 

12 

21.5585 

44 

13 

22.7974 

45 

14 

24.0070 

46 

15 

25.1882 

47 

16 

26.3418 

48 

17 

27.4684 

49 

18 

28.5688 

50 

19 

29.6437 

51 

20 

30.6940 

52 

21 

31.7202 

53 

22 

32.7232 

54 

23 

33.7038 

55 

24 

34.6626 

56 

25 

35 . 6004 

57 

26 

36.5179 

58 

27 

37.4159 

59 

28 

38.2951 

60 

29 

39.1563 

61 

30 

40.0002 

62 

31 

40.8276 

63 

FRCM  CALIBRATION 


(S) 


41.6391 

42.4355 

43.2176 
43.9862 
44.7419 
45.4855 

46.2177 
46.9393 
47.6510 
48.3536 
49.0478 

49.7343 
50.4139 
51.0874 
51.7554 
52.4187 
53.0781 

53.7343 
54.3881 
55.0401 
55.6911 
56.3419 
56.9932 
57.6457 
58.3003 
58.9575 
59.6183 
60.2832 
60.9531 
61.6287 
62.3108 
63.0000 


123 


APPENDIX  V 


Polynomial  Approximations  For 


and  C5^2(x) 


The  functions  are  defined  by  this  integral: 


C  (x)  =  ^7  f 

P  p '  J 

r  r  o 


2  2 
e  +x 


exp(-e)de 


The  integral  can  be  approximated  by  a  polynomial  with  a  high  degree  of  accuracy, 
for  the  required  range  of  x  as  shown  by  Burke  and  Hara  (9)*  For  the  integral 


c3/2(x) » 


C3/2(x) 


4  3  2 

x  +a^x  +^2X  +aix+a0 

x^+b  x^+b, x^+box^+b0x2+b  x+b 
5  4  3  2  1  0 


where 


ao 

= 

2.3983474 

X 

10 

a. 

= 

1.1287513 

X 

10 

1 

a2 

= 

1.1394160 

X 

io2 

- 

2.4653115 

X 

10 

3 

-2 

b0 

= 

1.8064128 

X 

10 

b 

- 

9.3877372 

1 

bn 

_ 

1.4921254 

X 

102 

2 

2 

2.8958085 

X 

10 

3 

b4 

= 

1.2049512 

X 

102 

N 

= 

2.4656819 

X 

10 
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For  the  integral  C^^2(x), 


'5/2 


(x)  = 


3  2 

x  +a^x  +ajX+a^ 


x5+b^x^+bgx3+b2x2+b1x+b0 


where 


aQ  =  1.1630641 
ax  =  1.6901002  x  10 
a2  =  6.6945939 
bQ  =  4.3605732 
bj  =  6.4093464  x  10 
b2  =  6.8920505  x  10 
b3  =  3.5355257  x  10 
b^  =  6.6314497 
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